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Solar  air  heaters  form  the  major  component  of  solar  energy  utilization  system  which  absorbs  the 
incoming  solar  radiation,  converting  it  into  thermal  energy  at  the  absorbing  surface,  and  transferring  the 
energy  to  a  fluid  flowing  through  the  collector.  The  efficiency  of  flat  plate  solar  air  heater  has  been  found 
to  be  low  because  of  low  convective  heat  transfer  coefficient  between  absorber  plate  and  the  flowing  air 
which  increases  the  absorber  plate  temperature,  leading  to  higher  heat  losses  to  the  environment 
resulting  in  low  thermal  efficiency  of  such  collectors.  Artificial  roughness  in  the  form  of  repeated  ribs  is 
the  most  effective  and  economic  way  of  improving  the  thermal  performance  of  solar  air  heater. 
This  paper  presents  an  extensive  review  on  the  research  carried  out  on  artificial  roughened  solar  air 
heater  ducts.  The  objective  of  this  paper  is  to  review  various  studies,  carried  out  on  thermal  as  well  as 
hydraulic  performance  of  artificial  roughened  solar  air  heater  ducts.  The  review  presented  in  this  paper 
will  be  useful  for  the  researchers  working  in  this  area. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Energy  is  a  basic  ingredient  needed  to  sustain  life  and  devel- 
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Nomenclature 

Re 

Reynolds  number 

Ra 

Rayleigh  number 

Ap 

Area  of  heated  plate,  m2 

Tf 

Average  temperature  of  air,  K 

cP 

Specific  heat  of  fluid  at  constant  pressure,  J / ( kg  K) 

Ts 

Temperature  of  the  sun,  K 

G 

Mass  velocity  of  fluid  through  the  collector,  kg/(m2  s) 

Ti 

Initial  temperature  of  air,  K 

Q 

Gap  or  discrete  distance,  m 

T0 

Final  temperature  of  air,  K 

Gp 

Gap  or  discrete  position,  m 

TP 

Average  plate  temperature,  I< 

Gr 

Grashof  number 

Ub 

Bottom  loss  coefficient,  W/(m2  1<) 

Gp/Lv 

Relative  discrete  or  gap  position 

Ue 

Edge  loss  coefficient,  W/(m2  I<) 

D 

Hydraulic  diameter  of  duct,  m 

UL 

Overall  heat  loss  coefficient,  W/(m2  K) 

e 

Rib  height,  m 

Ut 

Top  loss  coefficient,  W/(m2  K) 

e\D 

Relative  roughness  height 

V 

Velocity  of  air,  m/s 

fs 

Friction  factor  of  smooth  duct 

W 

Width  of  duct,  m 

f 

Friction  factor  of  roughened  duct 

W 

Width  of  V-rib,  m 

F 

Plate  efficiency  number 

W/w 

Roughness  width  ratio 

F0 

Heat  removal  factor 

AT 

Temperature  difference 

g 

Gap  or  discrete  width,  m 

AT/I 

Performance  parameter,  K  m2/W 

g/e 

Relative  gap  or  discrete  width 

H 

Depth  of  duct,  m 

Greek  letter  symbols 

h 

Convective  heat  transfer  coefficient,  W/(m2  I<) 

I 

Solar  intensity,  W/m2 

a 

Angle  of  attack,  degrees 

k 

Thermal  conductivity  of  air,  W/(m  K) 

n 

Collector  tilt  angle,  degrees 

L 

Length  of  test  section,  m 

>i 

Thermo-hydraulic  performance  parameter 

4 

Length  of  single  v-rib,  m 

p 

Density  of  air,  kg/m3 

m 

Mass  flow  rate,  kg/s 

Pm 

Density  of  manometric  fluid,  kg/m3 

Nu 

Nusselt  number  of  roughened  duct 

p 

Dynamic  viscosity  of  air,  N  s/m2 

Nus 

Nusselt  number  of  smooth  duct 

TOC 

Absorbance-transmittance  product 

P 

Pitch  of  the  rib,  m 

rjth 

Thermal  efficiency 

P/e 

Relative  roughness  pitch 

leff 

Effective  efficiency 

Pr 

Prandtl  number 

Qu 

Useful  heat  gain  rate,  W 

related  to  the  prosperity  or  the  standard  of  living.  Two  types  of 
energy  resources  are  available:  conventional  and  non- 
conventional  [1,2],  Conventional  energy  resources  such  as  fossil 
fuels  (coal,  crude  oil  and  natural  gas)  are  limited  in  amount.  Total 
energy  in  recoverable  conventional  energy  resources  is  estimated 
to  be  around  30-35 Q  (lQ=1018kJ)  while  the  global  energy 
consumption  rate  is  roughly  0.4-0.5  QJ yr.  Hence  conventional 
energy  resources  are  roughly  estimated  to  last  for  75-85  years. 
This  awareness  of  the  limited  nature  of  conventional  energy 
resources  gave  rise  to  the  search  of  alternate  energy  resources. 
Non-conventional  or  alternate  energy  resources  can  be  divided 
into  two  groups,  namely,  renewable  and  non-renewable  resources. 
Renewable  resources  are  those  which  have  a  short  period  of 
renewal  (upto  a  few  years)  such  as  solar  energy,  wind  energy, 
biomass  energy,  hydroenergy,  ocean  and  tidal  energy.  Solar  energy 
is  the  most  promising  of  all  these  alternatives  [3,4],  Solar  energy 
has  the  greatest  potential  among  all  the  sources  of  renewable 
energy  and  even  a  small  amount  of  this  renewable  source  of 
energy  is  sufficient  to  meet  the  total  energy  demand  of  the  world. 
If  we  can  use  5%  of  this  energy,  it  will  be  50  times  what  the  world 
will  require.  Solar  energy  is  readily  available,  well  distributed  and 
inexhaustible  for  all  practical  purposes,  and  has  no  polluting 
effects  upon  the  environment  when  converted  and  utilized.  Solar 
air  heaters,  because  of  their  inherent  simplicity  are  cheap  and 
most  widely  used  collector  devices.  Solar  air  heaters  are  being 
used  for  many  applications  for  low  and  moderate  temperatures. 
Some  of  these  are  space  heating,  crop  drying,  drying  of  concrete 
and  solar  dryer  [1,2, 5, 6, 7], 

The  use  of  artificial  roughness  on  the  underside  of  the  heated 
plate  can  substantially  enhance  the  thermal  performance  of  the 


solar  air  heater  due  to  increase  in  convective  heat  transfer 
coefficient  from  the  plate  to  air.  Surface  roughness  is  one  of  the 
first  techniques  to  be  considered  as  a  means  of  augmenting  forced 
convection  heat  transfer.  In  order  to  attain  higher  convective  heat 
transfer  coefficient  it  is  desirable  that  the  flow  at  the  heat  transfer 
surface  should  be  turbulent.  However,  the  turbulence  created  in 
the  core  can  increase  the  fan  power  exorbitantly.  It  is  therefore, 
desirable  that  the  turbulence  must  be  created  only  veiy  close  to 
the  heat  transfer  surface,  i.e.  in  the  laminar  sub-layer  only,  where 
the  heat  exchange  takes  place.  However,  as  pointed  out  above,  it  is 
necessary  that  while  creating  turbulence  to  break  the  laminar  sub¬ 
layer,  the  core  flow  should  not  be  disturbed  so  as  to  avoid 
excessive  losses.  This  can  be  achieved  by  using  artificial  roughness 
with  roughness  height  being  such  that  it  does  not  project  into  the 
core  but  is  of  the  height  that  just  project  out  of  laminar  sublayer. 
Numbers  of  experimental  investigations  involving  roughness  ele¬ 
ments  of  different  shapes,  sizes  and  orientations  with  respect  to 
flow  direction  have  been  carried  out  in  order  to  obtain  an 
optimum  arrangement  of  roughness  geometry.  Hans  et  al„  [8] 
carried  out  a  review  of  roughness  geometry  in  solar  air  heater 
ducts.  They  discussed  different  roughness  geometries  used  in  solar 
air  heater  ducts  and  explained  the  concept  of  artificial  roughness, 
effects  of  various  roughness  parameters  on  the  flow  pattern  and 
also  briefly  discussed  and  reviewed  the  roughness  geometries 
used  in  solar  air  heater  ducts.  The  objective  of  this  paper  is  to 
review  various  studies,  in  which  different  artificial  roughness 
elements  are  used  to  enhance  the  heat  transfer  rate  with  little 
penalty  of  increase  in  friction  losses  and  also  discuss  the  thermo- 
hydraulic  (thermal  as  well  as  hydraulic)  performance  of  artificial 
roughened  solar  air  heater  ducts. 
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Table  1 

Artificial  roughness  geometries  and  parameters  investigated  by  various  investigators. 


Sr.  no. 

Investigators 

Rib  geometry 

Parameters 

1. 

Prasad  and  Saini  [10] 

Transverse  continuous  rib 

P/e,  e/D 

2. 

Sahu  and  Bhagoria  [14] 

Transverse  broken  rib 

P/e,  e/D,  a 

3. 

Gupta  et  al.[15] 

Inclined  continuous  rib 

P/e,  e/D,  a 

4. 

Aharwal  et  al.[16] 

Inclined  rib  with  gap 

P/e,  e/D,  a,  d/W,  g/e 

5. 

Momin  et  al.[17] 

V-shaped  rib 

P/e,  e/D,  a 

6. 

Singh  et  al.[37] 

Discrete  v-rib 

d/W,  g/e,  P/e,  e/D,  a 

7. 

Karwa  et  al.[20] 

Chamfered  rib 

P/e,  e/D,  cp 

8. 

Saini  and  Saini  [21] 

Expanded  metal  mesh 

e/D,  S/e,  L/e,  a 

9. 

Saini  and  Saini  [22] 

Arc  shaped  rib 

e/D,  a 

10. 

Karmare  and  Tikekar  [23] 

Metal  grit  rib 

P/e,  e/D,  Ijs 

11. 

Bhagoria  et  al.  [24] 

Wedge  shaped  rib 

P/e,  e/D,  cp 

12. 

Saini  and  Verma  [25] 

Dimpled  shaped  rib 

P/e,  e/D 

13. 

Lanjeware  et  al.  [26] 

W-shaped  rib 

P/e,  e/D,  a 

14. 

Kumar  and  Bhagoria  [27] 

Discrete  w-shaped  rib 

P/e,  e/D,  a 

15. 

Varun  et  al.  [28] 

Combination  of  inclined  and  transverse  rib 

P/e,  e/D 

16. 

Jaurker  et  al.  [30] 

Rib-groove 

P/e,  e/D,  g/p 

17. 

Bopche  and  Tandale  [31] 

U-shaped  rib 

P/e,  e/D,  a 

18. 

Hans  et  al.  [32] 

Multi  v-shaped  rib 

P/e,  e/D,  a ,  W/w 

19. 

Patil  et  al.  [56-581 

Discrete  v-rib  with  combined  staggered  rib  pieces 

P/P,  S' IS,  r/e,  e/D.  P/e,  a 

20. 

Sethi  et  al.  [59] 

Dimpled  shaped  rib 

P/e,  e/D,  a 

21. 

Kumar  et  al.  [60,61] 

Multi  v-shaped  rib  with  gap 

G JLv ,  g/e,  W/w,  e/D,  P/e,  a 

2.  Artificially  roughened  solar  air  heater  ducts 

The  application  of  artificial  roughness  in  the  form  of  fine  wires 
or  ribs  having  different  geometries  on  the  heat  transfer  surface  has 
been  recommended  to  increase  the  heat  transfer  coefficient  by 
several  investigators.  The  use  of  artificial  roughness  on  the  under¬ 
side  of  the  absorber  plate  disturbs  the  viscous  sub-layer  of  the 
flowing  medium.  It  is  well  known  that  in  a  turbulent  flow  a  sub¬ 
layer  exists  in  the  flow  in  addition  to  the  turbulent  core.  The 
purpose  of  the  artificial  roughness  is  to  make  the  flow  turbulent 
adjacent  to  the  wall  in  the  sub-layer  region.  Artificial  roughness 
can  be  produced  by  machining,  casting,  forming,  welding  and 
other  methods. 

Artificial  roughness  can  be  produced  on  a  surface  by 

(i)  Blasting  sand  over  it. 

(ii)  Fixing  grooves  and  ridges. 

(iii)  Fixing  ribs  of  different  geometry  such  as  round,  rectangular  or 
V-shaped  or  Inclined  broken  ribs,  V-shaped  discrete  ribs  etc. 

(iv)  By  fixing  wires. 

(v)  By  fixing  three  dimensional  roughness  geometries. 

The  studies  on  the  effect  of  roughness  geometry  parameters  on 
the  heat  transfer  and  friction  factor  have  been  carried  out  to 
determine  the  most  efficient  configuration  and  optimum  dimen¬ 
sions  of  the  roughness  geometry  elements.  Some  of  the  important 
geometries  and  the  related  parameters  that  characterize  the 
roughness  geometry  and  substantially  influence  the  performance 
are  given  in  Table  1 . 

3.  Types  of  artificial  roughness  elements  investigated 

Prasad  and  Mullick  [9]  studied  the  effect  of  protruding  wires  on 
friction  factor,  heat  transfer  rate  and  plate  efficiency  factor  of  a 
solar  air  heater  used  for  drying  of  agricultural  products.  High  mass 
flow  rates  were  used  to  have  turbulent  flow  in  the  ducts  and 
protruding  wires  of  1  mm  diameter  tripped  the  laminar  sub-layer. 
They  compared  the  heat  transfer,  friction  factor  and  plate  effi¬ 
ciency  factor  of  roughened  corrugated  and  roughened  plane 
absorber  plates  with  that  of  corrugated  without  roughness  and 
smooth  absorber  plates.  The  protruding  wires  enhanced  the  heat 


transfer  rate  for  the  roughened  air  heaters  and  plate  efficiency 
factor  improved  from  0.63  to  0.72  resulting  in  14%  improvement  in 
the  performance. 

Prasad  and  Saini  [10,11]  studied  the  effect  of  roughness  and 
flow  parameters  such  as  relative  roughness  height,  e/D,  and 
relative  roughness  pitch,  P/e,  on  heat  transfer  coefficient  and 
friction  factor  for  a  fully  turbulent  flow.  They  found  that  the 
maximum  heat  transfer  occurred  in  the  vicinity  of  reattachment 
points  and  reattachment  of  free  shear  layer  does  not  occur  if 
relative  roughness  pitch,  P/e,  is  less  than  about  8-10.  Optimal 
thermohydraulic  (thermal  as  well  as  hydraulic)  performance  is 
achieved  for  roughness  height  slightly  higher  than  the  transition 
sub  layer  thickness.  For  relative  roughness  height,  e/D,  value  of 
0.033  and  relative  roughness  pitch,  P/e,  value  of  10,  maximum 
enhancement  in  heat  transfer  rate  and  friction  factor  was  reported 
to  be  2.38  and  4.25  times  respectively  in  comparison  with 
smooth  duct. 

Verma  and  Prasad  [12]  carried  out  an  outdoor  experimental 
investigation  for  thermohydraulic  (thermal  as  well  as  hydraulic) 
optimization  of  the  roughness  geometry  and  flow  parameters  for 
Reynolds  number,  Re,  range  of  5000-20,000,  relative  roughness 
pitch,  P/e,  range  of  10-40  and  relative  roughness  height,  e/D, 
range  of  0.01-0.03.  The  optimal  value  of  roughness  Reynolds 
number,  e+,  was  found  to  be  24  and  corresponding  to  this  value, 
thermohydraulic  efficiency  parameter  value  was  found  to  be  71 
percent.  Heat  transfer  enhancement  factor  was  found  to  vary 
between  1.25  and  2.08  for  the  range  of  parameters  investigated. 
Karwa  [13]  experimentally  investigated  the  effect  of  repeated 
rectangular  cross-section  ribs  on  heat  transfer  coefficient  and 
friction  factor  for  duct  aspect  ratio,  W/H,  range  of  7.19-7.75, 
relative  roughness  pitch,  P/e,  value  10,  relative  roughness  height, 
e/D,  range  of  0.0467-0.050,  Reynolds  number,  Re,  range  of  2800- 
15000.  It  was  explained  that  vortices  originating  from  the  rough¬ 
ness  elements  beyond  the  laminar  sub  layer  are  responsible  for 
heat  removal  as  well  as  increase  in  pressure  drop.  The  enhance¬ 
ment  in  the  Stanton  number  was  reported  to  be  65-90%  while 
friction  factor  was  found  to  be  2.68-2.94  times  over  smooth  duct. 

Sahu  and  Bhagoria  [14]  experimentally  investigated  the  effect 
of  90  broken  ribs  on  the  enhancement  of  the  thermal  perfor¬ 
mance  of  solar  air  heaters  for  fixed  roughness  height,  e,  value  of 
1.5  mm,  duct  aspect  ratio,  W/H,  value  of  8,  pitch,  P,  in  the  range  of 
10-30  mm  and  Reynolds  number,  Re,  range  of  3000-12000. 


Table  2 

Correlations  developed  for  heat  transfer  and  friction  factor  for  different  roughness  geometries  used  in  solar  air  heaters. 


Investigators 


Prasad  and 
Saini  [10] 


Verma  and 
Prasad 

[12] 


Gupta  et  al. 

[15] 


Aharwal 
et  al.  [16] 


Type  of  Geometry  Correlation  heat  transfer  coefficient  Correlation  friction  factor 

roughness 


Small  dia. 
Transverse 
protrusion  wire 


St  =  (£)  /  [l  +  y^{4.5  Pr°57(e+)0'28  -  0.95©°  53}J 


Nur  =  0.0859(P/e)- 0  054(e/D)0072Re0723  for  e  <  24 


Oblique  wire 
roughness 


Inclined  rib  with 
gap 


/  =  2/[0.95  (P/e)0  53 +  2.4  ln(D/2e)-3.75] 


fr  =  0.025(P/2)-°0206(e/D)0021  Re~’ 25 


/  =  0.1991(e/D)0196(W/H)~0093Re~0165  x  [exp( -0.0993(1 -a/70”)2)] 


/  =  0.071(Re)-°,33(P/e),s3>!  x  exp[-0.44(ln(P/e))2](d/W)-°43 

xexp[-0.14(lnd/W)2](g/e)-°052  x  (a/60)067  x  exp[0.12  (lng/e)2](e/D)069 
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Table  2  ( continued ) 


Investigators  Type  of  Geometry 

roughness 


Correlation  heat  transfer  coefficient 


Correlation  friction  factor 


Momin  et  al.  V  shaped  wire 
[17]  roughness 


Singh  et  al.  V-shaped  with 
[37]  gap 


Nu  =  0.067(e/D)0424Re0888  x  (a/60°)-°077[exp{-0.782(ln(a/60o))2)] 


Nu  — 2.36  x  10~3Re09°(P/e)3-50 

x(a/60)-°023(d/W)-0043(g/e)-0014(e/D)°47 
xexp(  -0.84(ln(P/e)2))exp(  -0.72(ln(a/60)2)) 
x  exp(  -  0.05(ln(d/W  )2)  )exp(  —  0.1 5(ln(g/e)2)) 


/  =  6.266(e/D)°  565Re0425  x  (a/60°)-°°93[exp{-0.719(ln(a/60  ))2)] 


/  =  4.13  x  10-2Re-°126(P/e)274(a/60)-°034(d/W)-0  058 
x(g/e)0031  (e/D)°  70exp(  -  0.685(ln  (P/e)2)) 

xexp( -0.93(ln  (a/60)2))  x  exp(  -0.058(ln (d/W)2))exp(-0.21(ln (g/e)2)) 


Kumar  et  al.  Discrete 
[27]  W-shaped  ribs 


Hans  et  al.,  V-shaped 
[8]  discrete 


Hans  et  al.  Expanded  metal 
[8]  mesh 


Saini  and  Arc  shaped  rib 
Saini  [22] 


-A 


—If-  M-| 


Nu  =  0.105  x  (Re)0 


'  x  (g)0453  x  (isT0™  x  exp [-0.59  x  (In^"))2 


/  =  5.68  x  (Re)-040  x  (e/D)  0  59  x  0  081  x  exp[-0.57  x  (ln(|,“))2] 


Nu 


0.0054Re,29x(l),34(|),,,2(g)°27 
xexp[-0.376(ln  (l-(^))2] 


/  =  0.717Re-2991  x  (f)00636(?)°0712 

x(g)01,3(£)-00093  xexp[-1.26(ln(l-(i))2] 


Nu  =  4  X  10-4Re'  22(e/D)0625  x  (S/10e)222(L/10e)266 

x[exp{- 1.25  (ln(S/10e))2)]  x  [exp(- 0.824  In  (L/10e)2)] 


/  =  0.815  Re-0361(L/e)0266  x  (S/10e)-°19(10e/D)°591 


Nu  =  0.001047Re13,86(e/D)°3772  x  (a/90)-°1198 


/  =  0. 14408Re  -°  17103(e/D)° 1765  x  (a/90)01185 
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Karmare 

and 

Tikekar 

[23] 


Bhagoria 
et  al.  [24] 


Saini  and 
Verma 

[25] 


Kumar  et  al. 

[27] 


Layek  et  al. 
[29] 


Metal  grit  rib 


Wedge  rib 
roughness 


Dimple  shape  rib 


Discrete 
Wshaped  ribs 


Chamfered  rib- 
grooved 


Nu  —  2.4  X  1(T3  X  Re1-3  x  (e/D)042(l/s)-0146  x  (P/e)-027 


Nu  =  1.89  x  10-4(e/D)0426Re121(P/e)294 

x[exp(-0.71  ln(P/e))2J(0/KJ)  0018[exp{  I  50f Inf 0/ 1 0  ))2}J 


Nu  =  5.2  x  1(T4  Re127{(P/e)315) 

x  [exp(  -2.12 )( log  ( P / e > ) 2 J ( t? / O ) 0033  x  [exp(-1.3)(log(e/D))2] 


Nur  =  0.105  X  (Re)0  873  x  (e/Dh)0453 

x(a/60  )-°081  x  exp[  — .59  x  (ln(a/60  ))2] 


Nu  =  0.00225  X  Re0  92  ( e/D)° 52  (P/e)1 72  x  exp[-0.22(ln  0)2] 

xexp[— 0.46(  ln(P/e)}2](g/P)  121 0124  x  exp[-0.74(  ln(g/P))2] 


f  =  15.55  X  Re-0  263  X  (e/D)091(//s)-°27  x  (P/e)-°51 


/  =  12.44  x  (e/D)099  x  Re-°18  x  (P/e)-°52  x  (a/10°)049 


/  =  0.642  Re0423 (P/e)- 0465  x  [exp(0.054)(log(P/e))2](e/Dr°°214 
x[exp(0.84)(log(e/D))2] 


fr  =  5.68  X  (Re)-040  X  (e/Dh)059  x  (a/60")-°081 
xexp[— 0.579  x  (ln(a/60  ))2] 


/  =  0.00245  x  Re0124(e/D)0365(P/e)432(g/P)-1 124 

xexp[-1.09{  ln(P/e))2]  x  exp[-0.68((g/P))2]  x  exp[O.OO50] 
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Table  2  ( continued ) 


Investigators  Type  of  Geometry 

roughness 


Correlation  heat  transfer  coefficient 


Correlation  friction  factor 


9  -  CHAMFER  ANGLE 
«  1  RIB  HEIGHT 
P  :  PITCH 
w  *  RIB  WIDTH 
a  =  ANGLE  OF  ATTACK 


Jaurker  et  al.  Rib  groove- 
130]  compound 

roughness 


Hans  et  al.  Continuous  Multi 
[32]  v-rib 


Patil  et  al.  Discrete  v-rib 
[56-58]  with  Combined 
staggered  rib 
pieces 


w w 


W/w  =  5  W/w=  10 


Nu  =  0.00206  Re°  936(e/D)°  349(P/e)3  318 
xexp[-0.868(  ln(P/e))2](g/P)' 108 
xexp[2.486{  ln(g/P))2  +  1.405{  ln(g/P)}3] 


Nu  =  3.35  x  10~5Re°92(e/D)°77  x  (W/w)°43(a/90r°49 
xexp(0.117  In  (W/w)2) 
xexp(  — 0.61  In  (a/90)2 )(P/er8  54 
xcxp(  -2  0407  In  (P/e)2) 


Nu  =  0.0089  Re0  97 

xexp[ap+{2.42  ln(f  )  +  1.19}2  +  0.11/l  +  {2.5  ln(p'/p)  +  1.41)2J 
+0.14[  ln(r/e)]°  71 


/  =  0.00127  Re-° 199 (e/D)°  585 (P/e)7 19 
x(g/P,0645  x  exp[  - 1  854{  ln(P/e)}2] 
xexp[1.513{  ln(g/P))2+0.8662{  ln(g/P)}3] 


/  =  4.47  x  10-4Re-°3188(e/D)-°73 

x(W/w)-°22(a/90)-°-39exp(0.52  In  (a/90)2)  x  (P/e)89exp(2.133  In  (P/e)2) 


/  =  0.09Re-°18 

xexp[aiS  +  {3.181n(f)  +  1.56}2  +  0.08/l  +  {2.61  ln(p'/p)  +  1.40)2] 
+0.17{  In  (r/e)25) 
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Sethi  et  al.  Dimpled  shaped 

[59] 


Kumar  et  al.  Multi  v-shaped 
[61]  with  gap  rib 


ui  DopnrftbcittabeFbK 


it' 

t-l-t 


Nu  =  7.1  x  10-3  x  (e/D)0362  x  (P/e)"0047  x  (a/60)"00048 
x  exp[  -  0.7792(ln(a/60))2] 


Nu  =  8.532  x  1 0  - 3 (Re)0932  (g ) 0175  ) 0506 

r  in  / r-  \  —0.0348 

x[exp(-0.0753(Ln(£))2]  (£) 
x  |exp( - 0.0653 (Ln ) ) 2 j  (f )  “ 0 0708 
x[exp(-0.223(Ln(§))2] 
x  (hi)  00239  [exp(0.1153(Ln(^))2J 
x(f),,96[exp(-0.2805(Ln(£))2] 


/  =  4.869  x  10_1(e/D)°2663  x  Re~0223  x  (P/er0059  x  (a/60)°°042 
xexp[-0.4801(ln(a/60))2] 


f  =  3. 1934  (Re)-03151  (g)026s(ffi)01132 
x[exp,0.0974(LnS))2]  (£)°°6,° 

x  [exp(-0.1065(Ln(5i))2]  (f  )~0'769 
x[exp(-0.6349(Ln(§))2] 
x(il)01553  [exp(-0.1527(Ln(|i))2] 
x (e )  0 7947  [exp|0.1486(Ln(£))2] 
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Roughened  heated  plate  increased  the  heat  transfer  coefficient  by 
1.25-1.4  times  as  compared  to  smooth  rectangular  duct  under 
similar  operating  conditions.  Corresponding  to  roughness  pitch,  P, 
value  of  20  mm,  maximum  value  of  heat  transfer  was  obtained 
that  decreased  on  the  either  side  of  this  roughness  pitch,  P,  value. 
Based  on  the  experimental  investigation,  the  thermal  efficiency  of 
roughened  solar  air  heater  was  found  to  be  in  the  range  of 
51-83.5%  depending  upon  the  flow  conditions. 

Gupta  et  al.  [15]  studied  the  effect  of  relative  roughness  height 
(e/D),  inclination  of  rib  with  respect  to  flow  direction  (a)  and 
Reynolds  number,  Re,  on  the  thermohydraulic  performance  of  a 
roughened  solar  air  heater  for  transitionally  rough  flow  region 
(5  <  e+  <  70).  For  a  roughened  solar  air  heater,  maximum 
enhancement  in  Nusselt  number  and  friction  factor  was  reported 
to  be  of  the  order  1.8  and  2.7  times  respectively  corresponding  to 
angle  of  inclination  values  of  60°  and  70°,  respectively.  Best 
thermohydraulic  (thermal  as  well  as  hydraulic)  performance  was 
reported  for  relative  roughness  height  (e/D)  value  of  0.023  and 
Reynolds  number,  Re,  value  of  14,000. 

Aharwal  et  al.  [16]  experimentally  studied  the  effect  of  width 
and  position  of  gap  in  an  inclined  discrete  ribs  having  square 
cross-section  on  heat  transfer  and  friction  characteristics  of  a 
rectangular  duct.  The  duct  had  an  aspect  ratio,  W/H,  of  5.84, 
relative  roughness  pitch,  P/e,  of  10,  relative  roughness  height,  e/D, 
of  0.0377,  angle  of  attack,  a,  of  60°,  relative  gap  width,  g/e,  range  of 
0.5-2  and  relative  gap  position,  d/W,  varied  from  0.1667  to  0.667 
for  Reynolds  number,  Re,  range  of  3000-18000.  For  the  discrete 
ribs  and  continuous  rib  roughened  ducts,  the  enhancement  in  heat 
transfer  was  reported  to  be  in  the  range  of  1.71-2.59  times  and 
1.48-2.26  times  respectively  under  similar  operating  conditions. 
The  maximum  values  of  heat  transfer  ratio  (Nu/Nus),  friction  factor 
ratio  (flfs)  and  thermohydraulic  parameter  obtained  were  corre¬ 
sponding  to  relative  gap  width  (g/e)  value  of  1.0  and  relative  gap 
position  (d/W)  value  of  0.25  for  the  range  of  parameters  investi¬ 
gated.  Particle  Image  Velocimetry  (PIV)  system  was  used  to 
visualize  the  effects  of  rib  geometry  on  the  flow  behaviour. 

Momin  et  al.  [17]  experimentally  investigated  the  effect  of  geome¬ 
trical  parameters  of  v-shaped  ribs  on  heat  transfer  and  friction  factor 
characteristics  of  a  rectangular  duct  of  a  solar  air  heater.  The 
investigation  covered  a  Reynolds  number,  Re,  range  of  2500-18000, 
relative  roughness  height,  e/D,  range  of  0.02-0.034  and  angle  of  attack, 
a,  range  of  30  degree-90  for  a  fixed  relative  roughness  pitch,  P/e,  of 
10.  Rate  of  increase  of  Nusselt  number  was  observed  to  be  lower  than 
the  rate  of  increase  of  friction  factor  with  an  increase  in  Reynolds 
number.  The  maximum  enhancement  of  heat  transfer  and  friction 
factor  as  a  result  of  providing  artificial  roughness  had  been  found  to  be 
2.30  and  2.83  times  respectively  over  the  smooth  duct  for  an  angle  of 
attack,  a,  value  of  60°.  It  was  reported  that  for  relative  roughness 
height,  e/D,  value  of  0.034  and  angle  of  attack,  a,  value  of  60°,  v-shaped 
ribs  enhanced  the  value  of  Nusselt  number  by  1.14  and  2.30  times  over 
inclined  ribs  and  smooth  absorber  plate  respectively. 

Varun  et  al.  [18]  carried  out  a  comparative  experimental  study 
of  augmented  heat  transfer  and  friction  in  a  rectangular  duct 
having  rectangular  cross-section  ribs  arranged  in  transverse, 
inclined,  v-shaped  continuous  and  v-shaped  discrete  pattern  for 
duct  aspect  ratio,  W/H,  range  of  7.19-7.75,  relative  roughness  pitch, 
P/e,  value  10,  relative  roughness  height,  e/D,  range  of  0.0467-0.050 
and  Reynolds  number,  Re,  range  of  2800-15000.  The  enhancement 
in  the  Stanton  number  over  the  smooth  duct  was  reported  to  be  in 
range  of  65-90%,  87-112%,  102-137%,  110-147%,  93-134%  and 
102-142%  for  transverse,  inclined,  v-shaped  up  continuous, 
v-shaped  down  continuous,  v-shaped  up  discrete  and  v-shaped 
down  discrete  rib  arrangement  respectively.  The  friction  factor 
ratios  corresponding  to  these  arrangements  were  found  as  2.68- 
2.94,  3.02-3.42,  3.40-3.92,  3.32-3.65,  2.35-2.47  and  2.46-2.58 
respectively 


Karwa  et  al.[19]  experimentally  studied  the  heat  transfer  and 
friction  in  a  high  aspect  ratio,  W/H,  rectangular  duct  with  repeated 
rectangular  cross-section  ribs  on  one  broad  wall  in  v-shaped 
discrete  and  discontinuous  patterns  at  angles  of  inclination  of 
45°  and  60°,  relative  roughness  pitch,  P/e,  value  of  10.63,  relative 
roughness  length,  B/S,  values  3  and  6  with  Reynolds  number,  Re, 
varied  from  2850  to  15500.  It  was  observed  that  corresponding  to 
relative  roughness  length,  B/S,  value  of  6,  the  highest  and  the 
lowest  values  of  Stanton  number  ratio,  St/Sts,  were  obtained  for 
60  v-shaped  down  discrete  and  60°  v-shaped  up  discontinuous 
ribs  respectively.  Friction  factor  ratio,  ///s,  values  were  found  to  be 
maximum  and  minimum  for  60  v-shaped  down  discrete  and  45° 
v-shaped  up  discrete  ribs  corresponding  to  relative  roughness 
length,  B/S,  value  of  3.  Stanton  number  and  friction  factor  values 
were  reported  to  be  higher  for  v-shaped  down  pattern  in  compar¬ 
ison  with  v-shaped  up  pattern  and  60°  v-shaped  ribs  performed 
better  than  the  45°  v-shaped  ribs.  Rib  pattern  had  a  strong  effect 
on  the  Stanton  number  and  friction  factor  with  discrete  ribs 
performing  better  than  other  rib  configurations  from  thermohy¬ 
draulic  performance  point  of  view. 

Karwa  et  al.  [20]  performed  an  experimental  investigation  of 
Nusselt  number  and  friction  factor  for  rectangular  ducts,  having 
aspect  ratio,  W/H,  in  the  range  of  4.8-12,  and  roughened  with 
repeated  integral  chamfered  ribs.  The  roughness  geometry  and  flow 
parameters  considered  for  the  investigation  were  Reynolds  number, 
Re,  range  of  3000-20000,  relative  roughness  height,  e/D,  range  of 
0.014-0.0328,  relative  roughness  pitch,  P/e,  range  of  4.5-8.5  and 
chamfer  angle,  0,  varying  from  - 15°  to  18°.  Stanton  number  and 
friction  factor  increased  with  increase  in  chamfer  angle  and  attained 
maximum  value  corresponding  to  chamfer  angle,  0,  and  value  of  15°. 
It  was  reported  that  Stanton  number  decreased  while  friction  factor 
increased  with  increase  in  aspect  ratio,  W/H.  As  compared  to  the 
smooth  solar  air  heater  duct,  the  presence  of  chamfered  ribs  on  one  of 
the  broad  wall  of  duct  yielded  up  to  about  two-fold  and  three-fold 
increase  in  the  Stanton  number  and  the  friction  factor  respectively  in 
the  range  of  parameters  investigated. 

Saini  and  Saini  [21]  carried  out  an  experimental  investigation  to 
study  the  effect  of  wire  mesh  roughened  absorber  plate  on  Nusselt 
number  augmentation  and  friction  factor  characteristics  of  solar  air 
heaters.  The  investigation  considered  relative  longway  length  of 
mesh,  Lfe,  in  range  of  25-71.87,  relative  shortway  length  of  mesh, 
S/e,  in  range  of  15.62-46.87,  relative  roughness  height,  e/D,  in  range 
of  0.012-0.039  and  Reynolds  number,  Re,  in  range  of  1900-13000.  It 
was  reported  that  the  maximum  heat  transfer  of  order  4  times  over 
the  smooth  duct  were  obtained  corresponding  to  angle  of  attack  of 
61.9  ;  relative  longway  length  of  mesh,  L/e,  value  of  46.87  and  relative 
shortway  length  of  mesh,  S/e,  value  of  25.  Maximum  value  of  friction 
factor  was  reported  for  angle  of  attack  of  72°,  relative  longway  length 
of  mesh,  L/e,  value  of  71.87  and  relative  shortway  length  of  mesh,  S/e, 
value  of  15.  Saini  and  Saini  [22]  studied  the  effect  of  arc  shaped  ribs 
on  the  Nusselt  number  and  friction  factor  of  rectangular  ducts  with 
Reynolds  number,  Re,  relative  roughness  height,  e/D,  and  relative  arc 
angle,  a,  varying  from  2000  to  17000,  0.0213  -  0.0422  and  0.3333- 
0.6666  respectively.  It  was  reported  that  relative  arc  angle,  a,  had  an 
opposite  effect  on  Nusselt  number  enhancement  and  friction  factor. 
With  decrease  in  relative  arc  angle,  a,  value,  Nusselt  number  value 
increased  while  friction  factor  value  decreased.  Enhancement  of 
Nusselt  number  and  friction  factor  was  reported  to  be  of  order 
3.6  and  1.75  times  respectively  over  smooth  solar  air  heater  duct  for 
relative  arc  angle,  a,  value  of  0.3333  and  relative  roughness  height, 
e/D,  value  of  0.0422. 

Karmare  and  Tikekar  [23]  experimentally  investigated  the  heat 
transfer  coefficient  and  friction  factor  characteristics  of  a  rectan¬ 
gular  duct  having  absorber  plate  roughened  with  defined  grid  of 
metal  ribs  of  circular  cross-section  and  the  investigated  roughness 
geometry.  The  investigation  considered  relative  roughness  height, 
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Fig.  2.  Effect  of  relative  roughness  height  (e/D)  on  thermal  performance  [35], 


e/D,  range  of  0.035-0.044,  relative  roughness  pitch,  P/e,  range  of 
12.5-36,  relative  grit  length,  Z/s,  range  of  1.72-1  and  Reynolds 
number,  Re,  range  of  4000-17,000.  Enhancement  in  heat  transfer 


Fig.  4.  Thermal  efficiency  versus  thermal  performance  for  different  value  of  Re  and 
e/D  [36], 

rate  was  found  to  be  187%  and  the  friction  factor  increased  by  213% 
and  optimum  performance  was  observed  for  relative  grit  length, 
1/s,  value  of  1.72,  relative  roughness  height,  e/D,  value  of  0.044,  and 
relative  roughness  pitch,  P/e,  value  of  17.5  for  the  range  of 
parameters  studied.  Bhagoria  et  al.  [24]  experimentally  studied 
heat  transfer  and  fluid  flow  characteristics  in  a  solar  air  heater 
with  absorber  plate  roughened  with  wedge  shaped  transverse 
integral  ribs.  The  investigation  encompassed  the  Reynolds  num¬ 
ber,  Re,  range  of  3000-18000,  relative  roughness  height,  e/D,  range 
of  0.015-0.033  and  rib  wedge  angle,  <£,  range  of  8-12°.  It  was 
reported  that  heat  transfer  rate  and  friction  factor  increased  by 
2.4  and  5.3  times  over  smooth  duct  in  the  range  of  parameters 
investigated. 

Saini  and  Verma  [25]  studied  the  effect  of  roughness  geometry 
and  operating  parameters  on  heat  transfer  coefficient  and  friction 
factor  in  a  roughened  duct  provided  with  dimple-shape  roughness 
geometry  for  the  range  of  Reynolds  number,  Re,  from  2000  to 
12,000,  relative  roughness  height,  e/D,  from  0.018  to  0.037  and 
relative  roughness  pitch,  P/e,  from  8  to  12.  For  the  range  of 
parameters  investigated,  Nusselt  number  was  found  to  be  max¬ 
imum  corresponding  to  relative  roughness  height,  e/D,  of  0.0379 
and  relative  roughness  pitch,  P/e,  of  10.  For  fixed  value  of  relative 
roughness  pitch,  P/e,  of  10,  friction  factor  attained  the  maximum 
and  minimum  values  corresponding  to  relative  roughness  height, 
e/D,  values  of  0.0289  and  0.0189  respectively.  Lanjewar  et  al.  [26] 
investigated  the  heat  transfer  and  friction  characteristics  of  artifi¬ 
cially  roughened  ducts  with  W-shaped  ribs.  The  range  of  para¬ 
meters  for  this  study  has  been  decided  on  the  basis  of  practical 
considerations  of  the  system  and  operating  conditions.  The  duct 
has  a  width  to  height  ratio  (W'/H)  of  8.0,  relative  roughness  pitch 
(P/e)  of  10,  relative  roughness  height  (e/Dh)  of  0.03375  and  angle  of 
attack  of  flow  (a)  of  30-75°.  The  heat  transfer  and  friction  factor 
results  have  been  compared  with  those  for  smooth  duct  under 
similar  flow  and  thermal  boundary  condition  and  thermo- 
hydraulic  performance  has  been  investigated. 

Kumar  and  Bhagoria  [27]  studied  heat  transfer  and  friction 
characteristics  of  a  rectangular  channel  having  discrete  w-shaped 
rib  roughness.  The  experiment  encompassed  the  Reynolds  num¬ 
ber,  Re,  range  from  3000-15000,  relative  roughness  height,  e/D, 
range  from  0.0168-0.0337,  relative  roughness  pitch,  P/e,  value  10, 
and  angle  of  attack,  a,  range  from  30°-75°.  It  was  reported  that 
maximum  enhancement  of  Nusselt  number  and  friction  factor  in 
comparison  to  smooth  duct  was  of  order  of  2.16  and  2.75  times 
respectively  for  an  angle  of  attack,  a,  value  of  60°.  Flow  separation 
and  secondary  flow  resulting  due  to  the  presence  of  discrete 
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Fig.  5.  Thermal  efficiency  as  a  function  of  Reynolds  number  and  relative  gap 
position  [37]. 


Fig.  6.  Thermal  efficiency  as  a  function  of  temperature  rise  parameter  and  relative 
gap  position  [37]. 


Fig.  7.  Thermal  efficiency  as  a  function  of  Reynolds  number  and  relative  gap  width  [37], 


w-shaped  rib  combined  to  provide  an  optimum  value  of  angle  of 
attack,  a. 

Varun  et  al.  [28]  carried  out  an  experimental  study  on  heat 
transfer  coefficient  and  friction  factor  characteristics  by  using  a 
combination  of  inclined  and  transverse  ribs  on  the  heated  plate  of 
a  solar  air  heater  with  Reynolds  number,  Re,  ranging  from  2000  to 
14  000,  relative  roughness  pitch,  P/e,  range  of  3-8,  relative 
roughness  height,  e/D,  value  of  0.030,  duct  aspect  ratio,  W/H, 


Fig.  8.  Thermal  efficiency  as  a  function  of  temperature  rise  parameter  and  relative 
gap  width  [37], 


value  of  10  and  roughness  height,  e,  value  of  1.60  mm.  For  relative 
roughness  pitch,  P/e,  value  of  8,  the  best  thermal  performance  was 
reported.  Layek  et  al.[29]  carried  out  an  experimental  investigation 
to  study  the  heat  transfer  coefficient  and  friction  factor  for 
repeated  transverse  compound  rib-groove  arrangement  on  heated 
plate  of  a  solar  air  heater.  Four  relative  rib-groove  positions,  g/P, 
values  of  0.3,  0.4,  0.5  and  0.6  were  investigated  for  fixed  relative 
roughness  height,  e/D,  and  relative  roughness  pitch,  P/e,  values  of 
0.030  and  10  respectively.  It  was  found  that  relative  roughness 
pitch,  P/e,  value  of  10,  relative  groove  position,  g/P,  value  of 
0.4  provided  about  2.42  and  2.6  times  increase  in  the  Nusselt 
number  and  friction  factor  respectively  for  entire  range  of 
Reynolds  number,  Re,  studied. 

Jaurker  et  al.  [23]  experimentally  investigated  the  heat  transfer 
and  fluid  flow  characteristics  of  rib-groove  roughened  rectangular 
duct.  The  effect  of  relative  roughness  pitch,  P/e,  relative  roughness 
height,  e/D  and  relative  groove  position,  g/P,  on  heat  transfer 
coefficient  and  friction  factor  has  been  studied.  The  presence  of 
rib-grooved  artificial  roughness  enhanced  Nusselt  number  and 
friction  factor  up  to  2.70  and  3.60  times  respectively  in  comparison 
to  smooth  absorber  plate.  The  maximum  heat  transfer  rate 
occurred  for  a  relative  roughness  pitch,  P/e,  of  about  6  and  relative 
groove,  g/P,  value  of  0.4.  Bopche  and  Tandale  [31]  carried  out  an 
experimental  investigation  to  study  the  heat  transfer  coefficient 
and  friction  factor  of  a  rectangular  duct  roughened  artificially  with 
U-shaped  turbulator.  The  experimental  encompassed  Reynolds 
number,  Re,  range  from  3800  to  18000,  ratio  of  turbulator  height 
to  hydraulic  mean  diameter,  e/D,  is  varied  from  0.0186  to  0.03986, 
relative  roughness  pitch,  P/e,  range  of  6.67  to  57.14  and  value  of 
angle  of  attack,  a,  is  kept  constant  at  90  degree.  The  U-shaped 
roughness  geometry  shows  appreciable  heat  transfer  enhance¬ 
ment  even  when  Reynolds  number  value  is  less  than  5000  where 
ribs  are  generally  inefficient.  At  Reynolds  number  value  of  3800, 
maximum  enhancement  in  Nusselt  number  and  friction  factor  are 
reported  to  be  2.388  and  2.50  respectively. 

Hans  et  al.  [32]  investigated  heat  transfer  coefficient  and 
friction  factor  characteristics  of  Multiple  v-shaped  ribs  roughened 
rectangular  duct.  The  investigation  encompassed  Reynolds  num¬ 
ber,  Re,  from  2000  to  20000,  relative  roughness  height  (e/D)  values 
of  0.019-0.043,  relative  roughness  pitch  (P/e)  range  of  6-12, angle 
of  attack  (a)  range  of  30  -75°  and  relative  roughness  width  (W/w) 
range  1-10.  Based  on  the  experimental  results,  correlations 
for  heat  transfer  coefficient  and  friction  factor  were  developed. 
A  maximum  enhancement  of  heat  transfer  rate  and  friction  factor 
due  to  presence  of  such  an  artificial  roughness  has  been  found  to 
be  6  and  5  times,  respectively,  in  comparison  to  the  smooth  duct 
for  the  range  of  parameters  considered.  The  maximum  heat 
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Table  3 

Relative  performances  of  rough  tubes  [39]. 


Case 

Description 

Constraints 

Parameter 

Relative  mass 

velocity  (G*) 

1 

Reduced  surface 

P 

A 

(f/fs )05 

^  (f/fs)05 

area 

Ps 

Qs 

As 

(St/Sts)15 

(St/Sts)05 

2 

Increased  heat 

p 

-A  -1 

K 

(St/Sts) 

C*  =  (f//s),/3 

transfer 

Ps 

As 

J<s 

</(f/L ) 

3 

Reduced  pumping 

Q 

-A  -1 

p 

(f/fs) 

C*  =  (St/Sts) 

power 

Us 

As 

Ps 

(/(St/Sts) 

FUX  IO-3 


Fig.  9.  Nusselt  number  ratio  versus  Reynolds  number  [42], 


transfer  enhancement  has  been  found  to  occur  for  a  relative 
roughness  width  (W/w)  value  6  while  friction  factor  attains 
maximum  value  for  relative  roughness  width  (W/w)  value  of  10. 

Kumar  et  al.  [60,61]  experimentally  investigated  Nusselt  num¬ 
ber  and  friction  factor  characteristics  of  Multiple  v-shaped  with 
gap  ribs  roughened  rectangular  duct.  The  investigation  encom¬ 
passed  Reynolds  number,  Re,  from  2000  to  20000,  relative  rough¬ 
ness  height  (e/D)  values  of  0.022-0.043,  relative  roughness  pitch 
(P/e)  range  of  8-12,  angle  of  attack  (a)  range  of  30°-75°,  relative 
gap  distance  (Cd/L„)  range  of  0.24-0.80,  relative  gap  width  (g/e) 
range  of  0.5-1. 5  and  relative  roughness  width  (W/w)  range  1-10. 


4.  Heat  transfer  and  friction  factor  correlations 


5.3.  Thermal  performance  representation 


The  widely  used  method  of  performance  presentation  given  by 
National  Bureau  of  Standard  (NBS)  and  the  American  Society  of 
Heating,  Refrigeration  and  Air  Conditioning  Engineering  [33], 

The  collector  efficiency,  r\ ,  is  a  measure  of  the  collector 
performance  and  is  defined  as  the  ratio  of  the  useful  heat  energy 
gain  over  a  time  period  to  the  incident  solar  radiation  over  the 
same  time  period. 


n  = 


IQudt 

ApfHRdt 


(1) 


where, 


H  is  the  rate  of  radiation  intensity  on  unit  area  of  surface  of  any 
orientation, 

R  is  the  factor  to  convert  radiation  to  that  on  the  plane  of  the 
collector. 


The  instantaneous  rate  of  useful  energy  collection  is  the 
difference  between  the  energy  due  to  incident  solar  radiation  that 
is  absorbed  by  the  absorber  plate  and  the  rate  at  which  the  energy 
is  lost  due  to  the  difference  of  the  plate  and  the  ambient 
temperatures.  The  useful  energy  collected  in  a  solar  air  heater 
can  be  expressed  mathematically  as 

Qu  =  Ap[I(ra)  —  Ul](Tp  —  Ta)]  (2) 

To  overcome  the  difficulty  of  determining  the  average  plate 
temperature,  TP,  a  parameter  called  plate  efficiency  or  collector 
efficiency  factor,  F',  is  introduced.  It  is  defined  as 

p  Actual  useful  energy  collected 

~  Useful  energy  collected  if  entire  collector  surface 
were  maintained  at  average  fluid  temperature 


So  it  would  be  convenient  to  express  Eq.  (2)  in  terms  of  inlet  or 
average  fluid  temperatures,  by  introducing  the  collector  efficiency 
factor,  F'as  follows: 

Qu  =  F’Ap[I(za)  -  Ui](Tf  -  Ta)]  (4) 


where 

T  (To  +  Tj) 

Tf  =  ^— 


(5) 


In  most  of  the  solar  air  heaters  roughened  artificially  for  heat 
transfer  enhancement,  correlations  have  been  developed  for  heat 
transfer  coefficient  and  friction  factor  in  terms  of  roughness 
geometry  parameters.  Some  of  the  correlations  developed  for  heat 
transfer  coefficient  and  friction  factor  are  given  in  Table  2. 


5.  Thermal  performance  of  solar  air  heater 

It  was  pointed  out  earlier  that  the  low  thermal  efficiency  of 
solar  air  heaters  can  be  improved  by  using  artificial  roughness  in 
the  form  of  different  shapes  fabricated  in  various  arrangements  to 
create  turbulence  near  the  wall  or  to  break  the  viscous  sub  layer. 
As  a  result,  increasing  the  heat  transfer  coefficient,  thermal 
efficiency  can  be  increased  but  at  the  same  time  creating  turbu¬ 
lence  requires  additional  energy  which  has  to  be  supplied  by  fan 
or  blower  at  the  expense  of  electrical  energy. 


For  collector  without  recycling,  Ta=Tt ,  the  Eq.  (4)  reduces  to 
Qu  =  F'Ap[/(ra)  -  U(]<T0-Ti)/2]  (6) 

The  rate  of  useful  energy  collected  is  also  expressed  by 
considering  enthalpy  rise  of  the  air  as 

Q_u  =  mCP(T0-Ti)  (7) 

Thermal  efficiency  can,  now  be  expressed  as 


Then  the  efficiency  of  the  flat  plate  collector  can  be  alterna¬ 
tively  expressed  as: 

Ijth  =  F'[(ra)  -  UL](T0-T,)/23] 


(9) 


Effective  efficiency,  %  Efficiency,  % 
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Fig.  13.  Effective  efficiency  as  a  function  of  Reynolds  number  and  relative  gap 
position  [37]. 


follows: 


»/th  =  Fo 


(za)  -  UL 


(12) 


Here  Fa  the  collector  heat  removal  factor  is  defined  as  the  ratio  of 
the  actual  heat  transfer  to  the  maximum  possible  heat  transfer 
rate  when  the  absorber  plate  is  at  the  inlet  fluid  temperature.  The 
quantity  Fa  is  equivalent  to  the  effectiveness  of  a  conventional  heat 
exchanger,  which  is  defined  as  the  ratio  of  the  actual  heat  transfer 
to  the  maximum  possible  heat  transfer.  The  maximum  possible 
useful  energy  gain  (heat  transfer)  in  a  solar  collector  occurs  when 
the  whole  collector  is  at  the  inlet  fluid  temperature;  heat  losses  to 
the  surroundings  are  then  at  minimum.  The  collector  heat  removal 
factor  times  this  maximum  possible  energy  gain  is  equal  to  the 
actual  useful  energy  gain  Qu. 


Qu  =F0Ap[I(ra)  —  Ui](Ti  —  Ta)] 


(13) 


Re 

Fig.  12.  Variation  of  effective  efficiency  at  different  Reynolds  number  and  Insola¬ 
tion  [36]. 


The  relation  between  the  collector  efficiency  factor  F  and  heat 
removal  factor  F0  is  given  as 


F  0  = 


mCp 


ApUi 


F'UlAp  ,  , 

exp<>~mcr^-1 


(14) 


In  particular  case  of  solar  air  heaters  without  recycling  and 
where  the  inlet  air  temperature  coincides  the  ambient  tempera¬ 
ture  (i.e.  T,-  =  Ta),  Eq.  (12)  reduces  to 

t]  =  F0(za)  (15) 


This  expression  of  efficiency  does  not  allow  the  real  operative 
parameters  to  be  shown. 

In  view  of  these  limitations,  Biondi  et  al.  [34]  proposed  the 
following  equation  for  efficiency  of  solar  air  heaters,  taking  in  air  at 
ambient  temperature  (Tj  =  Ta),  which  is  utilized  for  representing 
thermal  performance: 


'7th  =  fp[(T«)  —  LfL](T0  —  Ta)/2/] 


(10) 


1th  =  F  O 


(za)-UL 


(16) 


i1th  =  mCP(T0-Ti)/IAp] 


(ID 


where  F0  is  the  heat  removal  factor  related  to  outlet  temperature 
and  can  be  expressed  as 


The  thermal  efficiency  of  a  solar  collector  can  also  be  expressed 
in  terms  of  another  factor  called  the  heat  removal  factor,  F0,  as 


F  o  = 


mCp 

APUL 


exp 


f'UlAp  1 
mCP  f 


(17) 
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Temperature  rise  parameter,  K  m2/W 


Fig.  14.  Effective  efficiency  as  a  function  of  temperature  rise  parameter  and  relative 
gap  position  [37]. 


Fig.  16.  Effective  efficiency  as  a  function  of  temperature  rise  parameter  and  relative 
gap  width  [37]. 


Reynolds  number  (Re) 

Fig.  15.  Effective  efficiency  as  a  function  of  Reynolds  number  and  relative  gap 
width  [37], 

5.2.  Thermal  efficiency 

The  design  of  a  solar  air  heater  using  artificial  roughness  can  be 
based  on  the  best  thermal  performance  where  the  thermal 
performance  of  the  collector  is  expressed  in  terms  of  useful  energy 
gain  by  the  collector  per  unit  solar  energy  input  and  it  can  be 
determined  as  follows: 

os) 


5.3.  Enhancement  of  collector  thermal  performance 

As  studied  earlier,  the  use  of  artificial  roughness  on  the  under¬ 
side  of  absorber  plate  of  the  solar  air  heater  leads  to  considerable 
enhancement  in  the  heat  transfer.  These  result  in  similar  enhance¬ 
ment  in  the  thermal  efficiency  of  the  solar  collector  [35],  Fig.  1 
shows  that  there  is  substantial  improvement  in  the  thermal 
efficiency  of  solar  air  heater  and  that  this  enhancement  is  a  strong 
function  of  roughness  parameters. 

Saini  [35]  investigated  the  thermal  performance  of  solar  air 
heater  having  artificially  roughened  absorber  plate  by  fixing 
expanded  metal  matrix  on  it.  The  effect  of  roughness  and  operat¬ 
ing  parameters  has  been  examined  and  a  comparison  of  perfor¬ 
mance  of  roughened  solar  air  heater  with  that  of  conventional 
solar  air  heater  having  smooth  duct  was  made.  Thermal  perfor¬ 
mance  plots  for  roughened  as  well  as  corresponding  smooth  air 


Fig.  17.  Effective  efficiency  versus  Reynolds  number  for  roughened  and  smooth 
solar  air  heaters  [43], 


heaters  have  been  prepared  using  experimental  data  under  var¬ 
ious  operating  parameters. 

Fig.  2  shows  the  effect  of  relative  roughness  height,  e/D  on  the 
thermal  efficiency  of  solar  air  heaters.  It  indicates  that  with  an 
increase  in  relative  roughness  height,  the  thermal  efficiency 
increases.  Fig.  3  shows  the  effect  of  relative  longway  length  of 
the  mess,  L/e  on  thermal  efficiency.  These  figure  shows  that  the 
thermal  efficiency  increases  with  as  increase  in  mass  flow  rate  as 
also  with  decreasing  values  of  temperature  rise  parameter. 

Gupta  et  al.  [15]  experimentally  investigated  the  thermal 
efficiency  of  solar  air  heater  having  artificially  roughened  absorber 
plate  by  fixing  inclined  rib  on  it.  It  has  been  shown  that  the 
thermal  energy  gain  and  pumping  power  increase  with  increase  in 
flow  rate.  In  the  low  flow  rate  range,  thermal  energy  gain  increases 
sharply  with  increase  in  flow  rate  while  rate  of  pumping  power 
consumption  is  low,  resulting  in  higher  value  of  effective  effi¬ 
ciency.  Beyond  this  range  of  flow  rate,  the  rate  of  increase  of 
pumping  power  is  very  high  while  the  rate  of  useful  energy  gain 
becomes  nearly  constant  resulting  in  decrease  in  the  value  of 
effective  efficiency. 

Karmare  and  Tikekar  [36]  studied  the  thermal  efficiency  of 
solar  air  heater  having  artificially  roughened  absorber  plate  by 
fixing  metal  grit  ribs  on  it.  Fig.  4  shows  the  thermal  efficiency 
versus  (T0  —  T))//  for  different  values  of  Reynolds  number  and 
relative  roughness  height  of  roughened  solar  air  heaters.  It  is 
observed  that  for  a  particular  value  of  Reynolds  number  thermal 
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Fig.  18.  Entropy  versus  temperature  rise  parameter  for  P/e=12  and  e/D— 0.0498  [52], 


Fig.  19.  Entropy  generation  number  versus  temperature  rise  parameter  for  P/e=12 
and  e/D=0.0498.  AT//  [52], 


increased  from  0.20  to  0.35,  and  as  the  relative  gap  position  is 
increased  further  there  is  a  negligible  change  in  enhancement 
Figs.  7  and  8  show  the  variation  of  thermal  efficiency  with 
Reynolds  number  and  temperature  rise  parameter,  respectively,  at 
different  values  of  relative  gap  width.  It  has  been  concluded  that 
the  thermal  efficiency  is  the  highest  at  a  relative  gap  width  of 
1.0  and  decreases  on  both  sides  of  this  gap  width.  The  inter-rib 
local  heat  transfer  coefficient  on  the  downstream  side  of  the  rib- 
gap  is  influenced  by  the  flow  velocity  through  the  gap  and  area 
disturbed  due  to  gap  flow.  As  the  gap  width  increases,  the  flow 
velocity  reduces  and  the  area  disturbed  increases.  The  two  factors 
have  opposing  effects  on  local  heat  transfer  coefficient.  There 
exists  an  optimum  relative  gap  width  at  which  heat  transfer 
coefficient  and  hence  thermal  efficiency  is  the  highest.  For  the 
present  study,  the  thermal  efficiency  is  the  highest  corresponding 
to  a  relative  gap  width  of  1.0.  It  has  been  concluded  that  the 
enhancement  in  thermal  efficiency  is  much  more  as  the  relative 
gap  width  is  increased  from  0.50  to  1.0,  and  the  enhancement 
decreases  at  a  slower  rate  as  the  relative  gap  width  is  increased 
further. 


6.  Thermohydraulic  performance  of  roughened  solar  air 
heaters 

Studies  have  been  conducted  to  develop  methodologies  for 
assessing  the  net  heat  gain  in  energy  due  to  artificial  roughness  as 
augmentation  of  heat  transfer  vis-a-vis  enhancement  of  pumping 
power  requirement.  Efforts  have  been  made  to  optimize  such  a 
system  for  a  required  temperature  rise  if  the  heat  transfer  and  the 
friction  factor  data  are  available.  The  existence  of  optimum 
condition  for  a  given  rough  surface  was  proposed  by  Sheriff  and 
Gumley  [38],  It  was  shown  that  the  characteristics  of  a  roughened 
surface  necessary  to  give  the  minimum  pumping  power,  in  the 
channel  can  be  written  as 


efficiency  increases  with  increase  in  (T0-Ti)//  at  constant  relative 
roughness  height.  It  is  also  observed  that  thermal  efficiency 
increases  with  relative  roughness  height. 

Singh  et  al.  [37]  experimentally  investigated  the  thermal 
efficiency  of  solar  air  heater  having  artificially  roughened  absorber 
by  fixing  discrete  v-shaped  ribs  on  it.  The  variation  of  thermal 
efficiency  with  Reynolds  number  and  temperature  rise  parameter 
at  different  values  of  relative  gap  position  for  a  fixed  value  of 
insolation  and  other  roughness  parameters  is  shown  in  Figs.  5  and  6, 
respectively.  The  thermal  efficiency  increases  with  the  increase  in 
relative  gap  position  up  to  0.65  and  then  decreases  with  further 
increase  in  relative  gap  position.  The  inter-rib  local  Nusselt  number 
in  the  span-wise  direction  decreases  continuously  from  the  leading 
edges  to  the  V-shaped  apex  of  V-shaped  down  rib,  the  gap  causes  an 
increase  in  inter-rib  local  Nusselt  number  downstream  of  the  gap. 
Therefore,  depending  upon  the  position  of  the  gaps,  the  average  heat 
transfer  coefficient  varies.  If  the  gaps  are  created  toward  the 
V-shaped  apex  region,  it  will  help  in  increasing  the  local  heat  transfer 
coefficient  in  the  low  heat  transfer  coefficient  region;  therefore,  the 
average  heat  transfer  coefficient  and  hence  thermal  efficiency  is 
enhanced.  This  happens  until  the  relative  gap  position  of  0.65.  If  the 
gaps  are  placed  too  close  to  the  V-shaped  apex,  it  is  anticipated  that 
the  two  high  heat  transfer  regions  generated  on  the  downstream 
side  due  to  gaps  in  both  legs  of  V-shaped  rib  may  overlap,  and 
therefore  will  not  produce  a  similar  effect  of  enhancement  in  the 
average  heat  transfer  coefficient.  Thus,  placing  the  gaps  closer  to  the 
V-shaped  apex  (d/W)  0.65  results  in  a  lesser  increase  in  average  heat 
transfer  coefficient  and  hence  thermal  efficiency.  The  enhancement 
in  thermal  efficiency  increases  as  the  relative  gap  position  is 


(19) 


Performance  of  the  roughened  tubes  in  heat  exchanger  design 
relative  to  smooth  tubes  of  equal  diameter  was  studied  by  Webb 
and  Eckert  [39]  with  a  view  to  obtain  enhanced  heat  transfer 
capacity,  to  reduce  the  frictional  power  requirement  and  to  permit 
a  reduction  of  heat  transfer  surface  area.  They  stated  the  design 
objectives  for  heat  exchanger  having  tubes  with  inner  smooth 
surfaces  operating  at  fixed  total  rate  and  entering  fluid  tempera¬ 
ture.  The  design  objectives  are  summarized  as 

To  obtain  reduced  heat  transfer  area  for  equal  heat  duty  and 
pumping  power 

(Q/Qs  =  l  and P/Ps  =  l) 


To  obtain  increased  heat  capacity  for  equal  surface  area  and 
pumping  power 

(A /As  =  1  and  P/Ps=l). 


To  reduce  pumping  power  expenditure  for  equal  heat  transfer 
and  surface  area. 

(Q/Qs  =  l  and  A/As  =  1). 


The  heat  transfer  enhancement  QjQs  equals  to  the  relative  con¬ 
ductance  Jf/Ks=hA//tsAs  for  the  rough  and  smooth  tubes  heat 
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Temperature  rise  parameter,  (AT/I) 

Fig.  20.  The  effect  of  relative  roughness  height  on  augmentation  entropy  genera¬ 
tion  number  [53]. 


Temperature  rise  parameter.  (AT/I) 


Fig.  21.  The  effect  of  relative  roughness  pitch  on  augmentation  entropy  generation 
number  [53]. 

exchangers.  The  relative  conductance  can  be  written  as 

K/Ks  =  St/StsA/As.G/Gs  (20) 


The  corresponding  pumping  power  ratio  is 
P/Ps=f/fsA/As.(G/Gs)3 


(21) 


Webb  et  al.  [39]  specified  that  it  may  be  necessary  to  operate 
the  rough  tube  at  a  lower  velocity  than  smooth  tube  for  equal 
pumping  power.  If  the  heat  exchanger  mass  flow  rate  and  tube 
diameter  are  held  constant,  this  may  be  accomplished  by  increas¬ 
ing  the  number  of  rough  tubes  (N)  in  parallel  by  ratio, 

N/Ns  =  G/Gs. 


Eliminating  G/Gs  from  Eqs.  (20)  and  (21) 

K/Ks  _  St /Sts 
(P/PS)(A/AS)2'3  (///s)1/3 


Temperature  rise  parameter,  (AI/I) 

Fig.  22.  The  effect  of  relative  groove  position  on  augmentation  entropy  generation 
number  [53]. 


Temperature  rise  parameter.  (AT/I) 


Fig.  23.  The  effect  of  chamfer  angle  on  augmentation  entropy  generation  number  [53]. 


Fig.  24.  Variation  of  exergy  efficiency  with  Reynolds  number  for  various  roughness 
geometries  (at  e/d= 0.035,  /=1000  W/m2  and  UL=5  W/m2-K)  [54]. 


Eq.  (22)  provides  an  expression  containing  the  parameters  K/KSt 
PIPS,  and  A/As  defined  in  terms  of  the  ratios  of  friction  factor  and 
Stanton  number  ///s  and  St/Sts,  respectively. 


For  each  of  three  application  cases  stated  above,  two  of  these 
parameters  are  unity.  Table  3  given  the  results  of  Eqs.  (20)  and  (21 ) 
when  applied  to  each  of  the  three  heat  exchanger  applications. 
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Fig.  25.  Variation  of  exergy  efficiency  with  Reynolds  number  for  various  roughness 
geometries  (at  e\d= 0.035,  /=1000  W/m2  and  14=10  W/m2-I<)  [54]. 


Fig.  26.  Variation  of  exergy  efficiency  with  Reynolds  number  for  various  roughness 
geometries  (at  e/d= 0.035,  /=800  W/m2  and  (it=5W/m2-K)  [54]. 


Fig.  27.  Variation  of  exergy  efficiency  with  Reynolds  number  for  various  roughness 
geometries  (at  e/d= 0.035,  /=800  W/m2  and  l/L=10  W/m2-I<)  [54]. 


Fig.  28.  Exergy  efficiency  as  function  of  Reynolds  number  and  relative  gap  position 
(d/W)[  55], 


The  last  column  of  the  Table  3  gives  the  relative  mass  velocity 
G*=GS/G  calculated  from  Eqs.  (20)  and  (21)  required  to  satisfy  the 
constraints. 

Lewis  [40]  proposed  the  criteria  of  evaluating  the  thermo- 
hydraulic  performance  of  roughened  surfaces  by  defining  an 
optimization  parameter  called  thermo  hydraulic  performance 
parameter  q,  which  is  expressed  as  follows: 

V  =  (St/Sts)3 /(f/fs)  (23) 

The  roughness  that  yields  a  value  of  this  parameter  greater 
than  unity  is  useful  and  out  of  various  possible  sets,  the  one  that 
yields  the  maximum  value  of  this  parameter  should  be  used. 


6.1.  Effective  efficiency 


Cortes  and  Piacentini  [41  ]  defined  the  effective  efficiency  (/7eff) 
of  solar  air  heater  as  the  ratio  of  the  net  exergy  gain  to  the  energy 
available  as  insolation.  The  net  energy  gain  is  obtained  by 
subtracting  the  pumping  power  from  the  energy  gained  by  the  air 
flowing  through  the  collector.  Thus,  the  effective  efficiency  repre¬ 
sents  the  net  gain  per  unit  solar  energy  input.  It  is  expressed  as 


Q  u  —  (Pm/C) 

neS=  Mp 


(24) 


where  ‘C  (=^?;m?;tr(;th)  is  the  conversion  factor  accounting  for  net 
conversion  efficiency  from  primary  thermal  energy  to  mechanical 
energy  of  pumping  power. 


Temperature  rise  parameter,  K  m2/W 

Fig.  29.  Exergy  efficiency  as  function  of  temperature  rise  parameter  (AT/I)  and 
relative  gap  position  (d/W)  [55], 


t]f:  efficiency  of  the  fan, 

>/m:  efficiency  of  the  electric  motor, 

i]a:  efficiency  of  electrical  transmission  from  the  power 

plant,  and 

i/th:  efficiency  of  thermal  conversion  of  the  power  plant. 
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Fig.  30.  Exergy  efficiency  as  function  of  Reynolds  number  and  relative  gap  width 

igle)  [55], 


Temperature  rise  parameter,  K  m2/W 


Fig.  31.  Exergy  efficiency  as  function  of  temperature  rise  parameter  (AT/I)  and 
relative  gap  width  (g/e)  [55], 


The  value  of  ‘C  as  recommended  by  Cortes  and  Piacentini  [41  ] 
is  0.18  (typical  values  of  efficiency  factors  being:  ^=0.65,  t]m= 0.88, 
>/tr= 0.925  and  )7th=  0.344). 

Effective  efficiency  has  been  employed  by  several  investigators 
[37,42,43,51-53]  to  account  for  the  increase  in  mechanical  energy, 
higher  grade  energy,  expended  in  an  effort  to  increase  the  thermal 
performance  of  the  system  as  a  result  of  use  of  artificial  roughness. 

Karwa  et  al.  [42]  experimentally  investigated  the  thermal 
efficiency  of  solar  air  heater  having  artificially  roughened  absorber 
plate  by  fixing  chamfered  ribs  on  it.  The  enhancement  in  the 
Nusselt  number  due  to  the  artificial  roughness  is  found  to  be  about 
50-120%  as  shown  in  Fig.  9.  Again  the  greatest  enhancement  in  the 
Nusselt  number  has  been  observed  for  the  air  heater  with  the 
largest  relative  roughness  height.  The  enhancement  in  the  heat 
transfer  rate  has  been  found  to  result  in  a  lower  temperature  of 
the  roughened  heated  plate  and  in  the  increased  outlet  air 
temperature  as  compared  to  the  smooth  plate.  Thus,  the  heat  loss 
from  the  collector  to  the  surroundings,  which  is  proportional  to 
the  temperature  excess  of  the  heated  plate  over  the  ambient,  is 
lower  for  the  roughened  collector.  The  artificial  roughness  on  the 
heated  plate  that  results  in  the  desirable  increase  in  the  thermal 
efficiency  also  results  in  an  undesirable  increase  in  the  pressure 
drop  due  to  the  increased  friction.  The  friction  factors  for  the 
roughened  solar  air  heater  ducts  are  about  1.8-2.3  times  of  that  for 


the  smooth  duct  solar  air  heater  ducts  for  the  relative  roughness 
height  of  0.0197  and  about  3.1 -3.9  times  for  the  relative  roughness 
height  of  0.0441  as  can  be  seen  in  Fig.  10. 

Karmare  and  Tikekar  [36]  experimentally  investigated  the 
thermal  as  well  as  the  effective  efficiencies  of  solar  air  heater 
having  artificially  roughened  absorber  plate  by  fixing  metal  grit 
ribs  on  it.  Fig.  11  shows  the  rate  of  useful  energy  gain  and  the 
power  requirement  of  the  fan  as  a  function  of  Reynolds  number.  It 
is  observed  that  the  rate  of  increase  of  useful  energy  gain  is 
relatively  higher  at  low  range  of  Reynolds  number,  whereas  it  is  a 
bit  lower  at  higher  range  of  Reynolds  number.  But  the  rate  of 
increase  of  power  consumption  is  low  for  lower  range  of  Reynolds 
number  and  increases  relatively  at  high  rate  as  Reynolds  number, 
Re,  increases.  The  power  consumption  does  not  exceed  the  rate  of 
useful  energy  gain,  i.e.  the  net  energy  gain  rate  is  positive  and  it  is 
also  clear  that  at  higher  Reynolds  number,  the  rate  of  useful 
energy  collected  becomes  almost  constant  but  the  power  con¬ 
sumption  rises  steeply.  The  thermal  as  well  as  the  effective 
efficiencies  of  the  collector  as  a  function  of  Reynolds  number, 
Re,  are  shown  in  Fig.  12  which  shows  that  the  thermal  efficiency 
increases  with  Reynolds  number,  Re,  whereas  the  effective  effi¬ 
ciency  attains  a  maximum  value  at  a  particular  Reynolds  number 
and  thereafter  decreases.  This  shows  that  an  optimum  operating 
condition  does  exist  for  the  given  roughness  configuration  at 
which  the  effective  efficiency  is  maximum  for  a  particular  Rey¬ 
nolds  number. 

Singh  et  al.  [37]  experimentally  investigated  the  effective 
efficiency  of  solar  air  heater  having  artificially  roughened  absorber 
by  fixing  discrete  v-shaped  ribs  on  it.  Fig.  13  shows  that  that  the 
optimum  value  of  relative  gap  position  for  discrete  V-shaped  down 
rib-roughness  varies  with  variation  in  Reynolds  number.  For 
Re=22,000  the  effective  efficiency  is  the  highest  at  d/w= 0.65, 
and  for  22000  <  Re  >  27000  the  effective  efficiency  is  the  highest 
at  d/w= 0.20.  It  has  been  concluded  that  the  enhancement  in 
effective  efficiency  is  much  more  when  d/w  is  increased  from  0.20 
to  0.35  as  compared  to  when  d/w  is  increased  further,  beyond 
d/w= 0.65  the  enhancement  reduces.  It  is  also  observed  in  Figs.  13 
and  14  that  at  higher  Reynolds  number  or  at  lower,  AT//  the 
effective  efficiency  of  conventional  flat  plate  solar  air  heater  is 
more  than  the  roughened  solar  air  heaters.  The  friction  factor  ratio 
of  V-shaped  down  rib-roughened  ducts  to  that  of  smooth  ducts 
increases  with  the  increase  in  values  of  Reynolds  number.  Thus,  as 
the  Reynolds  number  increases,  the  pumping  power  for  rough¬ 
ened  duct  increases  at  a  much  higher  rate  than  that  of  smooth 
duct.  Beyond  a  certain  value  of  Reynolds  number,  the  pumping 
power  becomes  much  more  for  roughened  solar  air  heater.  Hence 
the  roughened  solar  heaters  have  low  effective  efficiency  at  high 
Reynolds  number.  A  similar  trend  is  observed  for  temperature  rise 
parameter,  as  shown  in  Fig.  14. 

Fig.  15  shows  the  plot  of  effective  efficiency  as  a  function  of 
Reynolds  number  for  a  different  relative  gap  widths.  It  is  seen  that 
for  Re  =  22000  the  optimum  relative  gap  width  is  1.0,  whereas  for 
Reynolds  number  ranging  from  22000  to  23000  and  from  23000 
to  28500  the  optimum  relative  gap  widths  are  1.5  and  2.0, 
respectively.  For  Reynolds  number  of  more  than  28500,  the 
smooth  conventional  solar  air  heater  has  the  highest  effective 
efficiency.  It  is  observed  from  Fig.  16  that  for  AT//= 0.0041  K  m2/W, 
the  smooth  conventional  solar  air  heater  has  the  highest  effective 
efficiency.  Beyond  AT//  of  0.0041  K  m2/W,  the  discrete  V-shaped 
down  rib-roughened  solar  air  heater  has  the  higher  efficiency.  For 
AT//  range  of  0.0041-0.0050  K  m2/W  and  0.0050-0.0055  K  m2/W, 
the  optimum  relative  gap  widths  are  2.0  and  1.5,  respectively. 
For  AT//=0.0055  K  m2/W  the  optimum  relative  gap  width  is  1.0. 
They  showed  that  the  considerable  enhancement  in  effective 
efficiency  as  g/e  is  increased  from  0.5  to  1.0,  and  as  the  g/e  is 
increased  further  the  enhancement  decreases  at  a  slower  rate. 


118 


A.  Kumar  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  37  (2014)  1 00-122 


Mittal  et  al.  [43]  carried  out  a  study  in  order  to  compare  the 
effective  efficiency  of  solar  air  heaters  having  different  roughness 
geometries.  It  is  seen  from  Fig.  17  that  among  all  the  roughness 
elements  investigated,  solar  air  heater  having  angled  ribs  of  low 
roughness  height  (e/D)  has  been  found  to  have  better  effective 
efficiency  in  higher  range  of  Reynolds  number  (more  than  12000). 
However,  in  lower  range  of  Reynolds  number  (less  than  12000)  is 
the  solar  air  heaters  having  expanded  metal  mesh  as  artificial 
roughness  element  have  better  efficiency.  Further  it  is  also 
observed  that  the  effective  efficiency  of  smooth  solar  air  heater 
is  better  than  the  roughened  solar  air  heaters  in  the  range  of  very 
high  Reynolds  number.  In  order  to  determine  the  uncertainty  in 
the  results  a  sensitivity  analysis  has  been  carried  out  as  shown  in 
Fig.  17.  In  the  lower  range  of  Reynolds  number  the  maximum  value 
of  uncertainty  was  found  to  be  as  2.5%  and  for  the  expanded  metal 
mesh  geometry  and  for  inclined  ribs  it  has  been  observed  as  0.8%. 
However,  in  the  higher  range  of  Reynolds  numbers,  the  maximum 
uncertainty  has  been  found  to  be  as  0.5%  and  0.4%  for  expanded 
metal  mesh  and  angled  ribs  respectively. 

6.2.  Exergetic  efficiency 

Altfeld  et  al.  [44]  proposed  a  method  for  represents  efficiency, 
based  on  the  second  law  of  thermodynamics  to  determine  the 
optimum  values  of  roughness  geometry  parameters.  All  energy 
flows  are  evaluated  in  terms  of  exergy.  Exergetic  efficiency,  defined 
as  a  ratio  of  net  exergy  flow  to  exergy  flow  associated  with  solar 
irradiation  on  collector  surface,  which  can  be  expressed  as: 

>1ex g  =  Ys  (25> 

Where  En  is  net  exergy  flow  rate  and  is  expressed  as 

E„  =  IApiim’h  ~  pm(  1  -  yc)  (26) 


1  is  the  global  irradiation  (W/m2) 

Pm  is  pumping  power 

rjc  is  the  Carnot  efficiency  [  =  (1  -Ta/T/)] 


Es  is  the  exergy  associated  with  solar  irradiation  on  collector 
surface  and  is  given  as 


Tsun  is  the  equivalent  temperature  of  the  Sun  as  a  blackbody 
( ss  5762  I<) 

Exergy  analysis,  based  on  second  law  of  thermodynamics,  of 
renewable  energy  systems  is  gaining  popularity  amongst  research¬ 
ers  and  scientists  as  it  deals  with  minimization  of  irreversibility 
and  can  indicate  the  possibilities  of  thermodynamic  improvement 
of  the  process  under  consideration.  Recently,  the  concept  of  exergy 
has  been  applied  by  many  scientists  [45-49,51]  to  optimize  the 
performance  of  solar  air  heaters. 

Suzuki  [45]  presented  an  exergy  balance  equation  for  a  solar 
flat  plate  collector  and  an  evacuated  tubular  collector.  It  was 
reported  that  both  types  of  solar  collectors  have  nearly  equal 
capabilities  in  exergy  gain  despite  large  differences  in  technologi¬ 
cal  efforts  and  expenses  to  produce  them.  Fujiwara  [46]  analyzed 
the  performance  of  solar  collectors  from  exergy  point  of  view  and 
stated  that  friction  process,  which  is  an  energy  conversion  from 
mechanical  energy  to  thermal  energy,  can  be  treated  as  exergy 
loss.  Said  and  Zubair  [47]  analyzed  the  performance  of  two  flat- 
plate  solar  collectors  and  a  photovoltaic  panel  based  on  first  as 
well  as  second  law  of  thermodynamics.  Second  law  efficiency 


values  of  17%  and  11%  were  obtained  for  photovoltaic  and  flat  plate 
collectors  respectively. 

Naphon  [48]  presented  a  mathematical  model  for  predicting 
the  heat  transfer  characteristics,  thermal  performance  and  entropy 
generation  of  a  double  pass  solar  air  heater  having  longitudinal 
fins.  Effect  of  height  and  number  of  fins  on  the  entropy  generation 
was  considered  and  it  was  found  that  with  increase  in  height  and 
number  of  fins,  entropy  generation  decreases  while  thermal 
efficiency  increases.  Kurtbas  and  Durmus  [49]  experimentally 
evaluated  the  energy  efficiency,  friction  factor  and  dimensionless 
exergy  loss  of  a  solar  air  heater  having  five  solar  sub-collectors  of 
same  length  and  width  arranged  in  series  in  a  common  case  for 
different  values  of  Reynolds  number.  It  was  reported  that  tem¬ 
perature  difference  (T0-T,)  of  air,  pressure  loss  and  collector 
efficiency  are  the  major  parameters  that  affect  the  exergy  loss. 
Ucar  and  Inalli  [50]  experimentally  compared  the  performance 
and  second  law  efficiency  of  flat  plate  solar  collector  having 
absorber  plates  of  different  shapes  and  orientations  and  reported 
that  the  least  value  of  second  law  efficiency  is  obtained  in  case  of  a 
conventional  collector.  Gupta  and  Kaushik  [51]  undertook  a 
comparative  study  to  evaluate  the  performance  of  solar  air  heaters 
roughened  with  different  roughness  geometries  on  the  basis  of 
energy,  effective  and  exergy  efficiencies.  It  was  found  that  artificial 
roughness  on  heated  surface  effectively  increased  the  efficiencies 
in  comparison  to  smooth  surface.  The  thermal  efficiency  in  general 
increased  in  the  following  sequence:  smooth  surface,  circular  ribs, 
v-shaped  ribs,  wedge  shaped  rib,  expanded  metal  mesh,  rib- 
grooved,  and  chamfered  rib-groove.  The  effective  efficiency  based 
criteria  also  follows  same  trend  of  variation  among  various 
considered  geometries,  and  trend  was  reversed  at  very  high  values 
of  Reynolds  number.  The  exergetic  efficiency  based  criteria  also 
follows  the  same  pattern;  however,  the  trend  was  reversed  at 
relatively  lower  value  of  Reynolds  number  and  for  higher  range  of 
Reynolds  number,  the  exergetic  efficiency  approaches  zero  or  may 
be  negative.  It  was  found  that  for  the  higher  range  of  Reynolds 
number,  circular  ribs  and  v-shaped  ribs  yield  higher  value  of 
exergy  efficiency  while  for  low  Reynolds  number  range,  chamfered 
rib-groove  yields  higher  value  of  exergetic  efficiency. 

Kumar  et  al.  [52]  experimentally  investigated  the  effective 
efficiency  of  solar  air  heater  having  artificially  roughened  absorber 
by  fixing  discrete  V-shaped  ribs  on  it.  Figs.  18  and  19  depict  the 
effect  of  relative  gap  position  on  the  entropy  and  entropy  genera¬ 
tion  number.  For  the  temperature  rise  parameter  lower  than 
0.038  I</m2W  relative  gap  position  of  0.15  has  minimum  value  of 
entropy  and  entropy  generation  number.  For  the  temperature  rise 
parameter  between  0.038  K/m2W  to  0.016  I</m2W  relative  gap 
position  of  0.35  has  minimum  entropy  and  entropy  generation 
number.  For  the  temperature  rise  parameter  higher  than  0.016  I</m2W 
relative  gap  position  of  0.15  has  minimum  value  of  entropy  and 
entropy  generation  number. 

Layek  et  al.  [53]  experimentally  investigated  the  second  law 
optimization  of  solar  air  heater  having  artificially  roughened 
absorber  by  fixing  chamfered  rib-groove  ribs  on  it.  Fig.  20  shows 
the  augmentation  entropy  generation  number,  Na  as  a  function  of 
temperature  rise  parameter  and  relative  roughness  height.  It  is 
evident  that  the  ability  of  the  chamfered  rib-groove  roughened 
absorber  plate  to  reduce  irreversibility  depends  on  the  thermo¬ 
fluid  operating  regime,  and  temperature  rise  parameter. 
The  entropy  generation  rate  decreases  sharply  with  increase  in 
temperature  rise  parameter  reaches  its  minima  for  temperature 
rise  parameter  of  about  0.007  I</m2W  and  increases  steadily  with 
further  increase  in  temperature  rise  parameter.  The  entropy 
generation  decreases  with  the  increase  in  relative  roughness 
height.  Fig.  21  depicts  the  effect  of  relative  roughness  pitch  on 
the  augmentation  entropy  generation  number,  Na.  For  the 
temperature  rise  parameter  higher  than  0.006  K/m2W  relative 
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Fig.  32.  Comparison  various  ribs  with  previous  investigations. 


roughness  pitches  of  6  have  the  best  performance  while  the 
relative  roughness  pitch  of  10  is  the  lowest.  For  the  temperature 
rise  parameter  lower  than  0.006  K/m2W  relative  roughness 
pitches  of  4.5  have  the  best  performance. 

Fig.  22  shows  the  effect  of  relative  groove  position  on  augmen¬ 
tation  entropy  generation  number,  Na.  It  can  be  seen  that  the 
relative  groove  position  of  0.4  has  the  best  performance  for  the 
temperature  rise  parameter  of  higher  than  0.004  I</m2W  and 
0.3  has  the  lowest.  For  the  temperature  rise  parameter  lower  than 
0.004  K/m2W  the  trend  is  reversed.  Fig.  23  shows  the  effect  of 
chamfering  of  the  rib  top  on  augmentation  entropy  generation 
number,  Na  with  variation  in  temperature  rise  parameter.  It  can  be 
seen  that  due  to  change  in  chamfer  angle  from  51  to  181  has  a 
remarkable  effect. 

Sharma  and  Varun  [54]  observed  the  variation  in  exergy 
efficiency  (r|Ex)  with  Reynolds  number  (Re)  for  different  values 
of  relative  roughness  height  (e/D).  It  is  seen  from  Figs.  24  to  27  that 
at  low  Reynolds  number  (less  than  12000),  the  exergy  efficiency 
increases  with  an  increase  the  values  of  relative  roughness  height 
(e/D).  However,  it  decreases  rapidly  at  higher  Reynolds  number 
(beyond  12000).  This  is  due  to  the  fact  that  the  rate  of  useful 
energy  collected  is  less  and  the  friction  losses  are  more  for  higher 
values  of  relative  roughness  height.  The  effect  of  irradiance  and 
overall  heat  loss  coefficient  are  found  to  be  is  similar  as  discussed 
earlier.  The  exergy  efficiency  (r|Ex)  also  follows  similar  trend,  for 
different  geometries,  as  shown  in  Figs.  3-6.  For  lower  range  of 
Reynolds  number  (up  to  8000),  rib-grooved  ribs  has  better 
exergetic  performance  as  compared  to  the  other  geometries.  For 
moderate  value  of  Reynolds  number  (up  to  12000),  solar  air  heater 
having  arc  shaped  wire  results  in  better  exergetic  efficiency. 
However,  beyond  12000  Reynolds  number  values,  smooth  air 
heater  is  found  to  be  better  than  that  of  a  roughened  solar  air 
heater. 

Singh  et  al.[55]  experimentally  compared  the  exergetic  effi¬ 
ciency  of  flat  plate  solar  collector  having  absorber  plate  of  different 
shapes  and  orientations  and  reported  that  the  optimum  value  of 
exergetic  efficiency  is  obtained  in  case  of  a  conventional  collector. 
It  can  be  seen  from  Fig.  28  that  the  highest  iju  is  obtained  at  d/w  of 
0.65  for  Re  <  13,000,  and  for  Re  >  15,300  the  smooth  conventional 
solar  air  heater  results  in  the  highest  i/n  For  Reynolds  number  in 
the  range  of  13,000  to  15,300  different  d/w  are  observed  to  result 
in  the  highest  r/u  depending  on  the  Reynolds  number.  From  Fig.  29 
it  is  observed  that  for  AT/I  <  0.008  I<m2/W  the  smooth  conventional 
solar  air  heater  results  in  the  highest  >j„  and  for  AT/I  >  0.0175  Km2/W 
the  ribs  having  d/w  of  0.65  results  in  highest  r/a.  For  AT/I  in  the  range 


Table  4 

Fixed  parameters. 


S.  no. 

Description 

Parameter 

Value 

1 

Collector  length,  m 

L 

1.5 

2 

Collector  width,  m 

W 

1.0 

3 

Thermal  conductivity  of  back 
insulated  sheet,  W/mK 

ki 

0.037 

4 

Thickness  of  back  insulated  sheet,  m 

U 

0.05 

5 

Density  of  air,  kg/m3 

P 

1.1541 

6 

Dynamic  viscosity  of  air,  kg/s-m 

P 

1.865  x  1CT5 

7 

Thermal  conductivity  of  air,  W/mK 

k 

0.02624 

8 

Transmittance-absorptance  product 

(ra) 

0.8 

9 

Emissivity  of  absorber  plate 

£P 

0.9 

10 

Emissivity  of  transparent  glass  sheet 

£g 

0.88 

11 

Ambient  temperature,  K 

Ta 

300 

12 

Wind  velocity,  m/s 

vw 

1.0 

of  0.008  to  0.0175  I<m2/W  the  highest  r|u  is  observed  for  different  d/w 
depending  on  value  of  AT/I.  The  enhancement  in  exergetic  efficiency 
is  maximum  at  d/w  of  0.65  except  at  low  AT/I.  The  enhancement  is 
more  when  d/w  is  increased  from  0.20  to  0.35  as  compared  to  when 
d/w  is  increased  further.  The  enhancement  reduces  at  a  slower  rate  as 
d/w  is  increased  further. 

Fig.  30  show  that  the  highest  i/u  is  obtained  at  g/e  of  1.0  for 
Re  <  12,500,  and  for  Re  >  15,200  the  smooth  conventional  solar  air 
heater  results  in  the  highest  r/n.  For  the  Reynolds  number  in  range 
of  12,500-15,200,  different  g/e  are  observed  to  result  in  the 
highest  i/ii  depending  on  the  Reynolds  number.  From  Fig.  31  it  is 
observed  that  for  AT/I  <  0.008  I<m2/W  the  smooth  conventional 
solar  air  heater  results  in  the  highest  t/u,  and  for  AT/I  >  0.0131  I<m2/W 
the  ribs  having  g/e  of  1.0  results  in  the  highest  t/u-  For  AT/I  in  the 
range  of  0.008  to  0.0131  I<m2/W  the  highest  t/u  is  observed  for 
different  g/e  depending  on  value  of  AT/I.  The  exergetic  efficiency  is 
maximum  at  g/e  of  1.0  except  at  low  AT/I.  There  is  considerable 
enhancement  in  the  exergetic  efficiency  as  g/e  is  increased  from0.5  to 
1.0  and  as  the  g/e  is  increased  further  the  enhancement  decreases  at 
a  slower  rate. 


7.  Comparison  effective  efficiency  of  roughened  solar  air 
heaters  duct 

The  solar  air  heater  having  artificially  roughened  absorber  plate 
have  higher  heat  transfer  coefficient,  through  accompanied  by  an 
increase  in  friction  factor.  Correlations  for  Nusselt  number  and 
friction  factor  have  been  developed  in  terms  of  roughness  and  flow 
parameters  by  various  investigators.  A  comparison  of  performance 
of  roughened  solar  air  heaters  with  that  of  conventional  solar  air 
heater  s  having  different  type  roughness  geometries  namely  here 
transverse  continuous  rib,  inclined  or  angled  rib,  discrete 
V-shaped  rib,  arc  shaped  rib,  metal  grit  rib,  V-shaped  rib,  single 
V-shaped  rib  with  staggered  rib  pieces,  wedge  shaped  rib, 
w-shaped  rib,  continuous  multi  V-shaped  rib,  discrete  multi 
V-shaped  rib  is  being  presented.  These  geometries  /shapes  have 
been  chosen  on  the  basis  of  the  fact  that  they  represent  relatively 
more  efficient  artificial  roughness  systems. 

7.2.  System  and  variable  parameters 

The  variable  parameters  comprise  of  roughness  geometry  i.e. 
relative  roughness  height  (e/D),  relative  roughness  pitch  (P/e), 
relative  width  ratio  (W/w),  relative  width  gap  (g/e)  etc.  and  design 
parameters  i.e.  temperature  rise  parameter  (AT /I)  and  average 
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solar  radiation  intensity  (I).  The  range  of  variable  parameters 
representing  roughness  geometry  has  been  selected  on  the  basis 
of  available  correlation  whereas  the  range  of  design  parameters  i.e. 
the  temperature  rise  parameter  and  insolation  has  been  decided 
on  the  basis  of  range  of  application  of  SAH  and  the  insolation 
available  at  a  place.  As  mentioned  in  Table  4  that  the  fixed 
parameters. 


7.2.  Procedure  for  prediction  of  effective  efficiency 

It  may  be  noted  that  since  the  objective  of  this  work  is  to 
present  a  methodology  for  optimal  design  of  the  solar  air  heater, 
the  results  need  to  be  presented  as  function  of  two  basic  design 
parameters  namely, 

•  Temperature  rise  parameter,  AT//  (ratio  of  air  temperature  rise, 
AT,  across  the  collector  to  the  average  intensity  of  solar 
insolation,  /) 

•  Solar  insolation  (/) 


The  calculation  starts  with  fixed  value  of  these  parameters 
(AT//  and  /)  and  proceeds  with  the  calculation  of  the  other 
parameters  for  a  given  collector.  Step  by  step  procedure  of 
calculation  is  given  below 

Step  1.  A  set  of  values  of  the  roughness  parameters  namely 
relative  roughness  height  (e/D),  relative  roughness  width 
(WD/WR),  wedge  angle  (<£),  relative  staggered  rib  posi¬ 
tion  (P/P),  relative  discrete  distance  (Gd//.„),  relative 
discrete  width  (g/e),  relative  staggered  rib  size  (r/e), 
relative  discrete  position  (S'/S),  angle  of  attack  (a), 
relative  roughness  pitch  (P/e)  and  angle  of  attack  (a)  is 
selected  for  which  the  calculation  is  to  be  performed. 

Step  2.  A  set  of  values  of  design  parameters  namely  temperature 
rise  parameter,  (AT//)  solar  insolation,  /  is  selected. 

Step  3.  Area  of  plate  is  calculated  as, 

Ap  =  WxL  (28) 


Step  4.  Hydraulic  diameter  of  duct  is  calculated  as; 
2  (WxH) 

W  +  H 


(29) 


Step  5.  Fixed  system  parameters  such  as; 

Thermal  conductivity  of  insulation,  /<, 

Thickness  of  insulation,  t, 

Transmittance-absorptance  product,  (m) 

Emissivity  of  the  absorber  plate,  ep 
Emissivity  of  the  glass  covers,  eg 
Fixed  operating  parameters  such  as; 

Atmospheric  air  velocity  (Vw)  and  atmospheric  tempera¬ 
ture  (Ta)  are  selected. 

Step  6.  The  temperature  rise  AT  of  air  and  outlet  temperature  T„ 
is  calculated  as, 

AT 

AT  =  — x  /  (30) 

T0  =  AT-t-T;  (31) 


Step  7.  Approximate  initial  mean  plate  temperature  is  assumed  as; 
T„+T, 

2 


TP  = 


-10°C 


Step  8.  Using  the  value  of  the  plate  temperature  Tp,  value  of 
top  loss  coefficient,  Ut  is  computed  by  using  equation 
given  as, 


of  Tp  +  Tg  )(TP  +  Tg) 

fkaNu\ 

((l/eP)  +  (l/eg)-l)  + 

V  Lg  )\ 

+  [(7£g(Tj,  +  Tg)(Tp  +  Tg)  +  hw]  1  +77-  (33) 

Kg 

where 


Tg 


(F,TP+cTa\ 

V  1+F1  ) 


where, 


[12  x  10  8(Ta  +  02TP)3  +  hw]  ^ +0.3fg 
r  1  — - 

[6  x  10-8(ep  +  0.028)(TP  +  0.5ro)3  +  0.6Ls  1 °2{(TP-T0)cos  /J)025] 

and 


f(Ts/Ta)  +  (hw/ 3.5)\ 
\  (l+(hw/3.5))  ) 


Ts  =  0.0522(Ta)1  5 


Nu  =  1  + 1.44[1  —  1708/Ra  cos  p]  + 

x{l  - 1708(  sin  1.8p)16/Ra  cos  p)+[(Ra  cos  p/5830)033  —  1] 

where,  Ra=Grx  Pr 

g/3'(TP  —  Tg)Lg 


[( To  +  T;)/2] 


Pr  = 


pCp 


l<a 


Back  loss  coefficient  Ub  is  expressed  as; 
Ub  =  kJU 


The  edge  loss  coefficient,  based  on  the  collector  area  AP  is 
given  as; 

,,  _  (L+W)teki 
e~  LWti 

finally, 

UL  =  Ut  +  Ub+Ue  (34) 


Step  9.  Useful  energy  gain  is  calculated  as; 
Qui  =  [Z(™)  —  Ui(Tp  —  Ta)]Ap 


(35) 


Step  10.  Mass  flow  rate  is  determined  from  the  expression  given 
as; 

™  Qui 


CpAT 


(36) 


Step  11.  Reynolds  number  of  flow  of  air  in  the  duct  is  computed  as; 

(37) 


D  GD 

Re  =  — 
P 


where,  G  is  the  mass  velocity  of  air  through  the  collector 
m 


G  = 


WH 


(38) 


2 


(32) 
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Step  12.  The  Nu  is  calculated  using  the  correlation  given  for  a 
particular  type  of  roughness  geometry.  The  correlations 
for  Nu  and  f  are  given  in  Table  3  for  the  roughness 
shapes  and  geometries  that  have  been  chosen  for  the 
present  investigation. 

Step  13.  The  convective  heat  transfer  coefficient  is  calculated  as 
follows, 


h 


Nu  k 
~ D~ 


(39) 


Step  14.  The  plate  efficiency  factor  is  then  determined  as, 


F'  = 


h 

h+lh 


(40) 


Step  15.  The  heat  removal  factor  is  calculated  as, 


_  mCp 
°  APUL 


F’UlM 
mCp  ) 


(41) 


Step  16.  The  useful  heat  gain,  QU2  is  computed  as, 

QU2  =  F0Ap[I(ra)  -  Ul(T0  -  Ti)]  (42) 


discussed  in  this  section.  The  results  of  the  roughened  solar  air 
heater  duct  are  also  compared  with  solar  air  heater  ducts  with 
smooth  surface  for  similar  AT//  values  and  numerical  conditions. 
Fig.  32  shows  the  effective  efficiency,  respectively,  as  functions  of 
AT//  for  various  roughness  geometries. 

Kumar  et  al.  [61]  showed  that  discrete  multi  V-ribs  on  the 
underside  of  the  absorber  plate  of  an  solar  air  heater  duct 
amplifies  thermal  as  well  as  hydraulic  efficiencies.  They  concluded 
that  the  higher  heat  transfer  rate  is  attributed  to  the  interaction  of 
secondary  flow  released  through  the  discrete  in  multiple  V-rib 
roughness  with  the  main  flow.  The  higher  thermal  performance  is 
attributed  to  the  interaction  of  the  secondary  flow  released 
through  the  discrete  with  the  main  flow.  This  increases  effective 
efficiency  through  the  discrete  width  area  behind  the  discrete. 
Discrete  multi  V-ribs  performed  better  than  non-discrete  multiple 
V-ribs.  As  shown  in  Fig.  36,  the  multi  V-with  gap  rib  geometry 
results  in  the  best  effective  efficienct  among  all  the  roughness 
shapes  investigated.  Tables  2  and  3  gives  the  range  of  optimum 
roughness  rib  parameters  and  thermohydraulic  performance  para¬ 
meters  compared  with  previous  studies. 


8.  Conclusions 


Step  17.  At  this  stage,  the  difference  between  the  two  values  of 
useful  heat  gain  QUl  and  Q„2  is  checked.  Ideally  the  two 
values  should  be  same.  However,  if  the  difference  in  two 
values  is  more  than  0.1%  of  QUl,  then  the  plate  tempera¬ 
ture  is  modified  as, 


Tp  —  Ta  4- 


7(ra)  — (QU2/Ap)’ 
UL 


(43) 


Step  18.  Using  the  new  plate  temperature,  the  calculations  from 
Step  6-17  are  performed  again  till  the  difference 
between  the  two  values  of  useful  heat  gain,  QU]  and  Q_U2 
is  reduced  to  a  value  below  0.1%  of  Qul. 

Step  19.  Friction  factor, /is  calculated  using  the  correlations  given 
for  a  particular  type  of  roughness  geometry.  The  correla¬ 
tions  used  for /are  given  in  Table  3  which  is  to  be  used 
for  analysis. 

Step  20.  Using  the  value  of  /  the  pressure  drop  (AP)d,  across  the 
duct  is  calculated  as  follow: 


(A  P)d 


4fLpV2 

2D 


(44) 


Step  21.  Using  the  value  of  (AP)d,  the  mechanical  power  (Pm)  is 
calculated  as, 


m(AP)d 


(45) 


Step  22.  Thermal  efficiency  is  calculated  as; 


'/th  =  F  o 


(t  a)- 


UL(T0-m 


I 


(46) 


Step  23.  The  effective  efficiency,  i/eff.  is  calculated  as; 


7 elf  — 


Q-u  —(F'm/C) 

lAp 


(47) 


The  effects  of  AT//  and  geometric  parameters  of  roughness  on 
effective  efficiency  in  an  solar  air  heater  duct  is  presented  and 


(i)  The  literature  review  shows  that  the  use  of  artificial  rough¬ 
ness  in  different  forms  and  shapes  is  an  effective  and 
economic  way  of  improving  the  performance  of  solar  air 
heaters. 

(ii)  Numbers  of  experimental  investigations  involving  roughness 
elements  of  different  shapes,  sizes  and  orientations  with 
respect  to  flow  direction  have  been  carried  out  in  order  to 
obtain  an  optimum  arrangement  of  roughness  element 
geometry. 

(iii)  Correlations  for  heat  transfer  and  friction  factor  have  been 
developed  which  are  applicable  to  wide  range  of  rib  config¬ 
urations  and  operating  parameters. 

(iv)  It  is  found  that  in  case  of  transverse  ribs,  two  fluid  vortices 
immediately  upstream  and  downstream  of  rib  are  stagnant 
on  account  of  which  fluid  temperature  as  well  as  wall 
temperature  increase  resulting  in  low  heat  transfer. 

(v)  Angling  of  transverse  rib  further  enhances  the  heat  transfer 
on  account  movement  of  vortices  along  the  rib  and  forma¬ 
tion  of  a  secondary  flow  cell  near  the  leading  end,  which 
results  in  local  wall  turbulence.  V-shaping  of  a  long  angled 
rib  helps  formation  of  secondary  flow  cells  as  compared 
single  transverse  ribs. 

(vi)  Creating  gap  in  V-shaped  rib  is  found  to  enhance  the  heat 
transfer  rate  by  breaking  the  secondary  flow  and  producing 
higher  level  of  turbulence  in  the  fluid  downstream  of  the  rib. 
Multi  v-shaped  rib  further  enhances  the  heat  transfer  rate 
because  increase  the  number  of  secondary  flow  cells 
several  times. 

(vii)  Methodology  of  artificial  roughness  and  experimental  stu¬ 
died  carried  out  by  various  investigators  have  been  discussed 
and  reported  in  detail.  Concept  of  thermal  efficiency,  ther- 
mohydraluic  performance  (effective  efficiency)  and  exergy 
efficiency  has  been  discussed.  It  is  observed  that  artificial 
roughness  is  a  good  technique  to  improve  the  thermo¬ 
hydraulic  performance  of  solar  air  heaters. 

(viii)  The  enhancements  in  the  Nusselt  number,  friction  factor  and 
thermal  efficiency  are  found  to  be  strong  functions  of  the 
relative  roughness  height.  The  greatest  enhancement  is 
observed  for  the  air  heater  with  the  highest  relative  rough¬ 
ness  height.  The  effective  efficiency  of  a  roughened  solar  air 
heater  increases  as  the  insloation  increases  for  Reynolds 
numbers  higher  than  10,000-12,000.  However  at  lower 
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Reynolds  number  the  thermohydraulic  efficiency  decreases 
with  increasing  insolation. 

(ix)  Most  of  the  technique  used  for  optimized  the  roughness 
geometries  used  in  solar  air  heater  for  optimum  values  of 
roughness  geometry  parameters  as  a  function  of  operating 
conditions  of  solar  air  heater  have  been  determined  on  the 
basis  of  thermal  efficiency,  effective  efficiency  and  exergetic 
efficiency  criteria. 
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